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This review covers latest developments in continuous-variable quantum-state tomography of opti- 
cal fields and photons, placing a special accent on its practical aspects and applications in quantum 
information technology. Optical homodyne tomography is reviewed as a method of reconstruct- 
ing the state of light in a given optical mode. A range of relevant practical topics are discussed, 
such as state-reconstruction algorithms (with emphasis on the majdmum-likelihood technique), the 
technology of time-domain homodyne detection, mode matching issues, and engineering of com- 
plex quantum states of light. The paper also surveys quantum-state tomography for the transverse 
spatial state (spatial mode) of the field in the special case of fields containing precisely one photon. 
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I. INTRODUCTION 
A. The concept of quantum tomography 

A quantum state is what one knows about a physical 
system. The known information is codified in a state 
vector or in a density operator p, in a way that 
enables the observer to make the best possible statis- 
tical predictions about any future interactions (including 
measurements) involving the syHtcni. Suc;h a definition 
has a comfortable interpretation within information the- 
ory, and so appears natural in the context of research in 
quantum information (QI). 

Imagine that an experimentalist, Alice, uses a well- 
characterized procedure to prepare an individual particle 
in a particular physical state. Since Alice possesses the 
information about the procedure, she can make definite 
predictions about the particle's behavior under various 
conditions, and is thus fully aware of the particle's state. 

Now suppose Alice sends the prepared particle to an- 
other party. Bob, who is not aware of the preparation 
procedure, but wishes to determine the state of the par- 
ticle. By making observations on the particle, Bob can 
obtain information about the physical state prepared 
by Alice by observing how it interacts with other well- 
characterized systems, such as a measurement appara- 
tus^. The amount and nature of this information de- 
pends strongly on whether the particle is macroscopic 
or microscopic. In the macroscopic, classical case. Bob 
can observe the individual particle's trajectory without 
disturbing it, and determine its state. 

In quantum mechanics, on the contrary, it is impossi- 
ble to learn the quantum state of any individual physical 
system. Each observation, no matter how subtle, will dis- 
turb its state just enough to prevent further observations 
from yielding enough information for a state determina- 
tion. This is the basis of quantum key distribution for 
cryptography (Bennett and Brassard, 1984). 

If Alice provides Bob with an ensemble of identically 
prepared systems, then he can measure the same variable 
for each system, and build up a histogram of outcomes, 
from which a probability density can be estimated. Ac- 
cording to the Born rule of quantum theory, this mea- 
sured probability density will equal the square-modulus 
of the state-vector coefficients, represented in the state- 
space basis corresponding to the measuring apparatus. 
This by itself will not yet yield the full state information, 
since the phase of the complex state-vector coefficients 
will be lost. 

As an example, measure the position x of each of 



^ We can interpret the quantum state as a belief, or confidence 

level, that a person has in his or her knowledge and ability to 
predict future outcomes concerning the physical system (Fuchs, 
2002). No measurements can, generally speaking, provide full 
information on Alice's preparation procedure. 



100,000 identically prepared electrons, which can move 
only in one dimension. This yields an estimate of the po- 
sition probability density, or the square- modulus |V'(a;)|^ 
of the Schrodinger wave function. If the wave function 
has the form |'!/'(a;)| exp[i(/)(a:)], where 4>{x) is a spatially 
dependent phase, then we will need more information 
than simply |V'(a;)|^ in order to know the wave function. 
If we are able to measure the momentum p of a second 
group of identically prepared electrons, then we can esti- 



mate the probability density 



where 



V'(p) = j V'( 



x) ex.p{—ixp/h)dx 



(1) 



is the Fourier transform of the spatial wave function. If 
we know a priori that the ensemble can be described 
by a pure state, then we can determine, by numerical 
methods, the complex wave function tp{x), up to certain 
symmetry transformations (such as a complex conjuga- 
tion) just from these two types of measurement. This is a 
classic example of phase retrieval (Gerchberg and Saxon, 
1972). 

In the typical case, however, we do not know ahead 
of time if the system's state is pure or mixed. Then we 
must make many sets of measurements on many sub- 
ensembles, every time modifying the apparatus so that 
sets of projection statistics associated with a different 
basis can be acquired. One can then combine these re- 
sults to reconstruct the density matrix of the state. The 
data do not yield the state directly, but rather indirectly 
through data analysis (i.e., a logical inference process). 
This is the basis of quanf,um- state tom,ography (QST). 
A set of observables whose measurements provide tomo- 
graphically complete information about a quantum sys- 
tem is called a quorum (Fano. 1957). 

Niels Bohr (1958) seems to have had an intuitive idea 
of QST, when he said, "A completeness of description 
like that aimed at in classical physics is provided by the 
possibility of taking cvny c;onc;eivable arrangement into 
account" . A more rigorous concept was developed in the- 
oretical proposals (Newton and Young, 1968; Band and 
Park, 1970, 1971, 1979; Bertrand and Bertrand, 1987; 
Vogel and Risken, 1989), followed by the first exper- 
iments determining the quantum state of a light field 
(Smithey et al., 1993a, b). Nowadays, quantum tomo- 
graphy has been applied to a variety of quantum systems 
and has become a standard tool in QI research (Paris and 
Rehacek, 2004). 

To continue the example of an electron moving in one 
dimension, a quorum of variables can be constructed by 
measuring different groups of electrons' positions x' after 
a variable length of time has passed. For example, in free 
space this is x' = x+pt/m, where m is the electron mass. 
The wave function at time t is 



j G{x',x;t)tp{x)dx, (2) 
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where 

is the quantum propagator appropriate to the wave equa- 
tion for the particle. The experimentally estimated prob- 
ability densities pr{x',t) = \tp{x',t)\^, for all t (positive 
and negative), provide sufficient information to invert 
Eq.(|2]) and determine the complex state function ^{x) 
(assuming the functions pr(x', t) are measured with very 
high signal-to-noise ratio). Note that Eq.(|2]) can be in- 
terpreted as a generalization of Eq.Q. As such, it cor- 
responds to a change of basis. 

If the state is not known beforehand to be pure (that 
is, the physical system's state is entangled with some 
other system), then it is described by a density matrix, 
p{xi, X2',t), in which case the probability densities corre- 
spond to the average, 

p{x'i,X2;t) = J J dxidx2 (4) 
X G*{x\,xi;t)G{x'2, X2;t)p{xi,x2\ 0). 

Through inversion of Eq. Q, the set of the measured 
probability functions pr(a;'; f) determines the density ma- 
trix p{xi,X2;Q)- 

This procedure works in principle for a Schrodinger 
equation with an arbitrary, known potential-energy func- 
tion. Such a method was proposed (Raymer et al, 
1994; Janicke and Wilkens, 1995; Leonhardt and Raymer, 
1996; Raymer, 1997b) and implemented (Kurtsiefer et 
al., 1997) for the transverse spatial state of an ensemble 
of helium atoms and the classical light beam (McAlister 
et al, 1995). 

B. Quantum tomography of light 

The current interest in QST is motivated by recent de- 
velopments in QI processing, which requires inter alia a 
technique for detailed characterization of quantum states 
involved (Paris and Rehacek, 2004). Additionally, signif- 
icant progress has been made in measurement technolo- 
gies, which now allow experimenters to measure a set 
of observables sufficiently diverse to allow reliable state 
reconstruction from the data. 

Among many physical systems in which QI processing 
can be implemented, light is of particular significance 
because it is mobile and thus irreplaceable as an infor- 
mation carrier in quantum communication networks. In 
this article, we review the basis of and methods for QST 
of optical fields. 

Even specialized to light, quantum tomography is too 
vast a field to be fully covered in a single review paper. 
Here we choose to concentrate on optical QST which in- 
volves measuring continuous degrees of freedom: field 
amplitude and/or spatial distribution. We study two 
tomographic problems which at first appear to have little 



in common. The first deals with the case in which the 
mode of the field is known (or chosen) a priori, and the 
state of this mode is to be determined. The second deals 
with the case that the full field is known to contain ex- 
actly one photon, but the manner in which this photon 
is distributed among spatial and spectral modes is to be 
determined. A careful analysis shows that these prob- 
lems are strongly related in their mathematical methods 
for reconstructing a quantum state. 



1. Optical homodyne tomography 

The first problem is the characterization of the state of 
the optical field in a certain spatiotemporal mode. The 
Hamiltonian of an electromagnetic mode is equivalent to 
that of the harmonic oscillator. Quantum states of light 
can thus be reconstructed similarly to motional states of 
massive particles discussed in the previous section. This 
is done by measuring quantum noise statistics of the field 
amplitudes at different optical phases (Leonhardt, 1997). 
The procedure of this reconstruction is known as optical 
homodyne tomography (OHT). 

It is interesting, in the context of this article, that ho- 
modyne tomography was the first experimental demon- 
stration of optical QST. Using balanced homodyne de- 
tection (BHD), Smithey et al. (1993a) measured a set 
of probability densities for the quadrature amplitudes 
of a squeezed state of light. These histograms were in- 
verted using the inverse Radon transform, familiar from 
medical tomographic imaging, to yield a reconstructed 
Wigner distribution and density matrix for a squeezed 
state of light. This 1993 paper introduced the term 
"tomography" into quantum optics. 

OHT is the subject of the next four sections of this 
article. In Sections II and III, we discuss the concept of 
homodyne tomography and the methods of reconstruct- 
ing the state's Wigner function and density matrix from 
a set of experimental data. Special attention is paid 
to the likelihood-maximization technique, which is now 
most commonly used. Section IV is devoted to technical 
issues arising in experimental OHT. In Section V, we dis- 
cuss applications of OHT in experiments of engineering 
and characterizing specific quantum states of light, such 
as photons, qubits, and "Schrodinger cat" states. 

In the context of applications, it is instructive to 
compare OHT with another technique: determining the 
quantum state of a system of dual-rail optical qubits^ 
by measuring relative photon number statistics in each 
mode and in their various linear superpositions. Due to 
its relative simplicity, this approach has been popular in 



^ In the dual-rail qubit, the logical value is assigned to a single 
photon being in one of two orthogonal modes A or B: = 

|1aiOs) , |l) = |0ai1b)i where the right-hand side is written 
in the photon number (Fock) basis for each mode. 
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a wide variety of experiments (see Altepeter et ai, 2004 
for a review). 

A textbook example of the above is the work of James 
et al. (2001). In this experiment the polarization state of 
a pair of entangled photons A and B generated in type-II 
parametric down-conversion was analyzed by measuring 
photon coincidence count statistics in sixteen polariza- 
tion projections. Tomographic analysis has revealed the 
photons to be almost perfectly in the state 

\^')^^{\HaVb) + \VaHb)), (5) 

where H and V indicate horizontal and vertical polariza- 
tion. It is tempting to say that such a state has high 
entanglement. 

This analysis does not reveal, however, that the pho- 
ton pair is generated not "on demand", but with some 
probability e^, which is usually low. A more complete 
representation of the state of the optical modes analyzed 
could be 



+ e{\\AH^Av'^BHQBv) 



(6) 



Pah'^av^br'^bv)) 



where, for example, indicates a one-photon state 

present in the vertical polarization mode of channel A. 
The bipartite entanglement is of the order |£^log(e)|, 
which is much less than 1. 

Eqs. ([5| and ^ reveal a significant limitation of the 
photon counting approach. This method works well if 
it is a priori known that the modes involved are in one 
of the qubit basis states or their linear combination. In 
practice, however, this is not always the case: photons 
can be lost, or multiple photons can be present where we 
expect only one. Such events compromise performance 
of quantum logical gates, but usually go unrecognized 
by the photon counting approach; they are simply elim- 
inated from the analysis'^. As a result, characterized is 
not the true quantum state of the carrier modes, but its 
projection onto the qubit subspace of the optical Hilbert 
space. This may lead to false estimation of gate perfor- 
mance benchmarks (van Enk et al., 2007). 

OHT, on the contrary, permits complete characteri- 
zation of the field state in a particular spatiotemporal 
mode, taking into account the entire Hilbert space of 
quantum optical states. It thus provides more reliable 
information about the performance and scalability of an 
optical QI processor. However, it is also more technically 
involved as it requires matched local oscillators, sophisti- 
cated detection electronics, larger POVMs and measure- 
ment data sets. It is thus less suitable for characterizing 



^ An important exception is the work by Chou et al. (2005), where 
the photon-counting approach is used, but the vacuum contri- 
bution is accounted for when evaluating the entanglement of a 
dual-rail qubit. 



multimode states: to date, the largest qubit systems mea- 
sured using photon counting contained six qubits (Lu et 
al, 2007) while with OHT, only single dual-rail qubits 
were reconstructed (Babichev et al., 2004b). 



2. Optical mode tomography 

An altogether different use of QST arises when a light 
field is known to contain a definite number of photons, 
but their distribution over spatial and/or spectral modes 
is unknown. If the set of modes is discrete, (e.g. in 
the case of polarization qubits), characterization can be 
done using the photon counting method discussed above 
(Altepeter et al., 2004). But if the distribution of light 
particles over electromagnetic modes is described by a 
continuous degree of freedom, methods of continuous- 
variable QST become irreplaceable. 

The problem of reconstructing the modal distribution 
of a field state is largely analogous to determining the 
spatial wave function V'(^ of a massive particle, as de- 
scribed by Eq.j2]). It turns out that the procedure out- 
lined via Eqs.( T|-([4]) also applies to QST for an ensem- 
ble of single photons, and is in fact quite similar to the 
method of OHT. This is consistent with the adoption 
of a sensible definition of a photon's spatial wave func- 
tion (Sipe, 1995; Bialynicki-Birula, 1996), in which the 
Schrodinger equation is replaced by the Maxwell equa- 
tions (since a photon cannot be strictly localized in space, 
some subtleties must be taken into account). In Section 
|VI| of this review we analyze various techniques and re- 
cent experimental progress in reconstructing spatial op- 
tical modes of single photons as well as entangled pairs. 



II. THE PRINCIPLES OF HOMODYNE 
TOMOGRAPHY 

A. Balanced homodyne detection 

The technique of balanced homodyne detection (BHD) 
and homodyne tomography has been extensively de- 
scribed in the literature, for example, in the textbook 
of Leonhardt (1997) and in recent reviews by Raymer 
and Beck (2004), Zhang (2004) and Zavatta et al. (2005). 
Here, we present only a brief introduction with concen- 
tration on theoretical aspects of mode matching between 
the local oscillator and the signal field. 

Figure [T] illustrates balanced homodyne detection, 
which is a means to measure the amplitude of any phase 
component of a light mode. In BHD, the weak signal field 

Es(t) (which may be multimode) and a strong coherent 

local oscillator (LO) field EL{t) are overlapped at a 50% 
reflecting beam splitter, and the two interfered fields are 
detected, temporally integrated, and subtracted. 

The signal electric field operator is written as a sum of 
positive- and negative-frequency parts, which are conju- 
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FIG. 1: Balanced homodyne detection 



.(+) 



gates of one another, Eg = Eg + Eg . The positive- 
frequency part can be decomposed into plane waves ac- 
cording to Dirac's quantization scheme 



Es ir,t) = ^ 



E 



flUJj 

2^ 



bjij exp(ikjf ~ iujjt), (7) 



where ujj, kj, and ej are, respectively, the mode fre- 
quency, wave vector, and the unit polarization vector; the 
creation and annihilation operators obey the commuta- 
tor [SjjSj,] — Sjji and are defined in some large volume 
V (which may be taken to infinity later) . It is convenient 
to consider the signal field in the paraxial approximation 
with z being the propagation axis. In this case ujj ~ ckjz 
and the polarization ej is along either x or y. 

The LO field is treated classically, and at each photo- 
diode face (z — 0) is assumed to be a strong coherent 
pulse propagating along the z axis. 



Ei+\r,t) 



2^ 



ULeLVLix, y) gL{t) exp{ik 



LZ 



(8) 
iujLt), 



where the coherent-state amplitude is = jaij e , and 
VL{x,y)g]^{t) is the normalized spatiotemporal mode. 

The local oscillator and the signal fields meet at a beam 
splitter, where they undergo the transformation 



(El^Es) 



El + Es El — Es 



V2 



(9) 



The difference of the numbers of photoelectrons recorded 
in the two beam splitter outputs is then, assuming a per- 
fect detection efllciency (see Raymer et ai, 1995; Raymer 
and Beck, 2004 for details and more general considera- 
tions) , 



SoV 



{2ELEs)dtdxdy 



cfiuj 

Det At 



(10) 



where the integration is done over the detector sensitive 
area and the measurement time At. Assuming that the 
above fully accommodate the local oscillator pulse, all in- 
tegration limits in Eq. (10 1 can be assumed infinite. The 
photon creation operator a) associated with the detected 
spatiotemporal mode is given by 



(11) 



where the C/s equal the Fourier coefficients for the LO 
pulse. 



L^3 



v*Lix,y)gl{t) 



(12) 



X exp (ikjxX -\- ikjyy — ic{kjz — kL)t) dt dx dy. 

When using a pulsed LO field EL{t), the concept of 
a light mode needs to be generalized beyond the com- 
mon conception as a monochromatic wave. As first dis- 
cussed by Titulaer and Glauber (1966), a polychromatic 
light wave packet can be considered a mode with a well 
defined spatial-temporal shape, whose quantum state is 
described in the usual way using photon creation and 
annihilation operators. For example, a one-photon wave- 
packet state is created by jla) = a} |vac) (more on this 



in Sec. VIA I. The meaning of Eqs.(lll, (12) is that the 



BHD detects the state of the electromagnetic field in the 
spatial-temporal mode defined by the LO pulse (Smithey 
et at, 1993a; Raymer et al, 1995; Raymer and Beck, 
2004). This allows temporal and spatial selectivity, or 
gating, of the signal field (not the signal intensity) . This 
gating technique (linear-optical sampling) has applica- 
tion in ultrafast signal characterization (Dorrer et al, 
2003; Raymer and Beck, 2004). 

As usual, the mode's annihilation operator can be 
expressed as a sum of Hermitian operators a = 
{Qe + iPe) /V^, called quadrature amplitudes, with^ 
[Qe,Pe] — i- For zero phase, Qg, Pg are denoted Q, P, 
respectively (so Qg — Q cos 9 + P sin 9) , and are analo- 
gous to position and momentum variables for a massive 
harmonic oscillator. For the LO phase equal to 9, BHD 
measures the quadrature amplitude 



(13) 



According to quantum mechanics, the probability density 
for observing the quadrature equal to Qg for the field in 
the signal mode given by the density operator p is 



piiQg,9)^{Qg\p\Qg) 



(14) 



where IQg) is the quadrature eigenstate. These proba- 
bility densities, also known as marginal distributions, are 



^ Some authors use the convention [Qg, Pg] = i/2. All quadrature- 
dependent plots used in this article have been (re-)scaled to com- 
ply with the uniform convention [Qg, Pg] = i. 
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histograms of the field amplitude noise samples measured 
with the homodyne detector. The optical phase plays the 
role of time in Eqs. ([iJ-Q. When it is varied over one 
complete cycle, quadrature amplitudes Qe form a quo- 
rum for QST (Vogel and Risken, 1989). 

In a practical experiment, the photodiodes in the ho- 
modyne detector are not 100% efficient, i.e. they do 
not transform every incident photon into a photoelec- 
tron. This leads to a distortion of the quadrature noise 
behavior which needs to be compensated for in the recon- 
structed state. We present, without derivation, a general- 
ization of the above expression for detectors with a non- 
unitary quantum efficiency rj (Raymer and Beck, 2004; 
Raymer et ai, 1995): 



\/27rcr2 



(15) 



where 2a^ ~ 1 /rj and the brackets indicate a quantum 
expectation value. The double dots indicate normal op- 
erator ordering (annihilation operators to the right of 
creation operators). 



B. Wigner function 

Because the optical state reconstructed using tomo- 
graphy is generally non-pure, its canonical representation 
is in the form of a density matrix, either in the quadrature 
basis or in the photon- number (Fock) basis. In the case 
of homodyne tomography it is convenient to represent 
the reconstructed state in the form of the phase-space 
quasiprobability density, the Wigner function (Wigner, 
1932). 



1 

2^ 



{Q+^Q'\p\Q-^Q')e-''''^ dQ'. 

(16) 

This object uniquely defines the state and, at the same 
time, is directly related to the quadrature histograms 
(14 15) measured experimentally (Raymer et at., 1995; 



Raymer and Beck, 2004) via the integral 

-\-OC +00 



(17) 



j j 5{Qg-Q cose -P sine) WDctiQ,P)dQdP 

oo —CO 
oo 

J Wuct {Qe COS e — Pg sine, Qe sin e + Pg cos e) dPg . 



In other words, the histogram pr((5g,0) is the integral 
projection of the Wigner function onto a vertical plane 
oriented at angle 6* to the Q axis (Fig. [2]). The "detected" 
Wigner function W^ot corresponds to the ideal Wigner 
function (16 1 for a loss- free detector, and for a detector 




FIG. 2: The Wigner function. The experimentally measured 
field quadrature probability density pT{Qg,6) is the integral 
projection of the Wigner function W{Q,P) onto a vertical 
plane defined by the phase of the local oscillator. 



via a convolution (Leonhardt and Paul, 1993; Kuhn et 
ai, 1994; Raymer et ai, 1995; Leonhardt, 1997) 



+ 00 +00 



WBet{Q,P) 



J W{Q',P') 



xexp 



7r(l - 77) 

— 00 —00 

(.Q-Q'Vv? + {P-P'Vv?' 



1-r? 



(18) 



dQ'dP'. 



III. RECONSTRUCTION ALGORITHMS 

A homodyne tomography experiment yields a set of 
pairs {Qrm ^m)i which can be binned up to form marginal 
distributions pr((3,6') for several local oscillator phases. 
Our next task is to develop mathematical methods that 
can be used to convert the experimental data into the 
state's density matrix and/or Wigner function. This is 
the subject of this section of the review. 

Mathematical methods of OHT can be divided into two 
categories. The so-called inverse linear transform tech- 



niques (Sec. Ill A I use the fact that the experimentally 



measured marginal distributions are integral projections 
of the Wigner function. Because integration is a linear 
operation, one can reverse it and reconstruct the Wigner 
function from the set of marginals in a procedure that 
somewhat resembles solving a system of linear equations 
of the form (17 1. We discuss this and related methods in 



Sec. Ill A they are reviewed in more detail in (Paris and 
Rehacek, 2004; Leonhardt, 1997; D'Ariano, 1997; Welsch 
et al, 1999; Raymer and Beck, 2004). 



For reasons discussed later (Sec. IIIB 



with quantum efficiency 77 it is obtained from the latter 



inverse linear 

transform methods are rarely used in modern OHT. More 
frequently, we employ methods of statistical inference, 
whose classical versions have been developed in tradi- 
tional statistics and data analysis (Paris and Rehacek, 
2004, chapters 2,3,6, and 10). A popular method is like- 
lihood maximization (MaxLik) , which looks for the most 
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probable density matrix that will generate the observed 
data. It is discussed in detail in Sec. lIIIBl Another statis- 
tical inference method, entropy maximization, is briefly 
reviewed in Sec. IIIICI 



State reconstruction via inverse linear 
transformation 



Inverse Radon transformation 



The projection integral ([17|), known as the Radon 
transform (Herman, 1980), can be inverted numerically 
using the back-projection algorithm, familiar from med- 
ical imaging (Herman, 1980; Leonhardt, 1997) to recon- 
struct the phase-space density WBct{Q,P)- 



1 



TT +00 



27r2 

-oo 

X K{Q cose + P sine -Qg)dQed9, 



(19) 



with the integration kernel 



K{x) 



+ 00 



(20) 



where V denotes a principle value integration. 

The kernel is infinite at a; = 0, so in numerical im- 
plementations of the inverse Radon transformation it is 
subjected to low pass filtering: the infinite integration 
limits in Eq. (20 1 are replaced by ifc^, with chosen 



so as to reduce the numerical artifacts associated with 
the reconstruction while keeping the main features of the 
Wigner function [see e.g. Fig. [sj^a)]. This method is 
known as the filtered back-projection algorithm. 

This strategy was used in the first QST experiments 
(Smithey et ai, 1993a; Dunn et al, 1995). In later im- 
plementations of this algorithm (Lvovsky and Babichev, 
2002), the intermediate step of binning the data and 
calculating individual marginal distributions associated 
with each phase was bypassed: the summation of Eq. ( 19 ) 
was applied directly to acquired pairs (0m,Qm): 



1 



N 



cos e^n + P sin e„i - Qm), 
(21) 

with phases uniformly spread over the 27r interval. 

A nonclassical state of light^ , after undergoing an op- 
tical loss, becomes nonpure. Therefore, typically, the re- 
constructed state is not pure, that is, Tr[/5Qj,j] ^ 1. A 



^ See, for example, Lvovsky and Shapiro (2002), as well as Zavatta 
et al. (2007) for a review of definitions and measurable criteria 
of a nonclassical nature of a state of light. 



special case is that of a coherent state, which is not non- 
classical and remains pure under losses. For such a state, 
one can reconstruct the Schroedinger wave function or a 
state vector, as demonstrated in Smithey et al. (1993b). 

A more general method for reconstructing the Wigner 
function of a state that has undergone optical losses is 
proposed in Butucea et al. (2005), where Eqs. (19 1 and 



( 20 1 are modified so as to incorporate the effect of nonuni- 
tary efficiency. This paper, as well as Gu^a and Artiles 
(2006), also perform a minimax analysis of the error in 
the evaluation of the Wigner function via the inverse 
Radon transformation. 

Given the experimentally reconstructed Wigner func- 
tion, we can inverse Fourier transform Eq. (16 1 to com- 



pute the density operator in the quadrature basis, and, 
subsequently, in any other basis. This scheme was ap- 
plied to reconstruct photon-number statistics {n\p\n), 
as well as quantum-phase statistics for squeezed and for 
coherent light (Beck et al., 1993; Smithey et al., 1993b; 
Smithey et al., 1993c). This calculation can however be 
significantly simplified, as discussed below. 



2. Pattern functions 

If the goal is to reconstruct the density operator of the 
ensemble, we can exploit the overlap formula 



+00 +00 



Tr(pi) ^ 2tt 



WpiQ,P)WAQ,P)dQdP (22) 



valid for any operator A and the associated density ma- 



trix W^{Q, P) as defined by Eq. ( 16 1 with p replaced by 

A. For example, given Amn — \'^){'>A with \m) and \n) 
being the Fock states, we write pmn = Tr(|m)(n| p) and 



use Eq. (22 1 to determine, one by one, the elements of 



the density matrix in the Fock basis. 

The intermediate step of reconstructing the Wigner 
function can however be sidestepped using an improved 
deterministic scheme introduced by D'Ariano et al. 
(1994), and refined several times to the present opti- 
mal form (Leonhardt et al., 1996, 1997; D'Ariano et al, 
2004). We combine Eqs. ([l9l and ([22]) to write 



TT +00 

Tr(pi) = / / w{Qe,0)F^{Qg,e)dQede (23) 

-oo 

where averaging is meant in the statistical sense over all 
acquired values of {Qe, e), and 



(24) 



+ 00 

- JJ K{Q COS e + Psm9~ Qe)Wj^{Q, P) dQ dP 
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is the sampling function. Given a specific operator A, 
the function does not depend on the experimental 
histogram pi{Qg,9), but only on the operator itself. It 
thus need be calculated only once, prior to the experi- 
ment, and substituted into Eq. ( 23 ) once the data become 
available. 



Specializing to the Fock basis, F„in{Q,6) = 
(l/7r)e*('"-")«M„„(Q), with Af™„(Q) being so-called 
pattern functions (D'Ariano et ai, 1994; Paul et ai, 
1995; Lconhardt and Raymer, 1996): 



+ 00 



-V 



1pTn{x)i'n{x) 

{Q-xY 



(25) 




where 



6 8 10 12 14 

n 



FIG. 3: Quantum optical state estimation from a set of 14152 
experimental homodyne measurements (Lvovsky and Mlynek, 
2002) by means of the inverse Radon transformation and the 
pattern-function method (a) and the likelihood maximiza- 
tion algorithm (b). The Wigner function and the diagonal 
elements of the reconstructed density matrix are shown. The 
inverse Radon transformation in (a) was performed by means 
of the filtered back-projection algorithm. The statistical un- 
certainties in (b) were determined by means of a Monte-Carlo 
simulation (see text). 



1/4 



Hn{x)_ 

\/2"n! 



exp 



X 



(26) 



are the Fock state wave functions — that is, wavefunc- 
tions of energy eigenstates of a harmonic oscillator. _ff„ 
denote the Hermite polynomials. Figure [3](a) shows an 
example of calculating the density matrix using the pat- 
tern function method. 

Efficient numerical algorithms for computing the pat- 
tern functions were given by Lconhardt et al. (1996) and 
Leonhardt (1997). In our experience, the most practical 
algorithm involves the irregular wave functions (/3„(a;), 
which are alternative, non-normalizable solutions of the 
time-independent Schrodinger equation for the harmonic 
oscillator. These functions obey a recursion 



^/2r^T2 



[xifr,{x) - ip'n{x)] (27) 



with 



^o(a:)-7r3/4exp(^--jerfi(a;) (28) 

and are all readily expressed through the error function 
erfi(x). Once the desired number of the irregular wave 
functions have been calculated, the pattern functions are 
obtained using 

/ N _ / d[^„,{x)Lpr,{x)]ldx for n>m , . 
M^n{x) - <^ a[V'„(a;)(^„(a;)]/9x for n< m. ^^^> 

The pattern function method can be extended to di- 
rect sampling, or quantum estimation (Paul et al., 1995; 
Munroe et al., 1995), of the expectation value of any 
operator directly without first reconstructing the state. 
In many cases, this requires fewer probability functions 
to be measured, since less complete information is being 
asked for. Indeed, since Wigner functions are linear with 
resp ect to their generating operators, we conclude from 
Eq. (24) that for any operator A — 'Y^ Amn 



FA{Qe,d) = ^F™„(Q,0)imn (30) 

E{n\A \m) Mmn{q) exp[i{m - n)9]. 



The expectation value (A) — Tr{pA) can then be calcu- 
lated according to Eq. (|23|. 



9 



For example, if we desire the photon-number probabil- 
ity prO), we choose Aj = \j) Then F^,{Q,9) — 
{l/Tr)Mjj{Q), which is independent of phase 9, and 
Eq.(|23| becomes (Munroe et al, 1995) 

pr(j) = Tr{p\j) 01) = / dQM,,{Q){w{Q.e))e (31) 

This is a convenient result, since only a single probability 
function need be measured, while sweeping or randomiz- 
ing the phase. Demonstrations of this technique can be 
found in Munroe et al. (1995), Schiller et al. (1996), as 
well as Raymer and Beck (2004). 

Although the techniques of OHT have been generalized 
to fields involving more than one optical mode (spatial, 
polarization, or temporal) (Raymer et al, 1996; Opatrny 
et al., 1996; D'Ariano et al., 2000), their practical appli- 
cation is challenging. This is a particular case where 
direct sampling is handy. We can apply it if our task 
is to determine the expectation values of certain observ- 
ables, but full reconstruction of the multimode state is 
not necessary. An example is the acquisition of correlated 
photon number statistics of two-mode fields (McAlister 
and Raymer, 1997b; Vasilyev et al., 2000; Blansett et al, 
2001; Voss et al, 2002; Blansett et al, 2005). 

B. Maximum-likelihood reconstruction 

1. Why maximum Ukelihood? 

Quantum state reconstruction can never be perfect, 
due to statistical and systematic uncertainties in the esti- 
mation of the measured statistical distributions. In both 
discrete- and continuous- variable domains, inverse linear 
transformation methods work well only when these un- 
certainties are negligible, i.e. in the limit of a very large 
number of data and very precise measurements. Other- 
wise the errors in the "right-hand sides" of the system 
of linear equations we are trying to solve can lead to in- 
accurate, even seemingly unphysical, features in the re- 
constructed state. For example, negative values may be 
found on the diagonal of the reconstructed density matrix 
and its trace is not guaranteed to equal one [Fig. [sj^a)]. 

In the case of continuous-variable tomography, there 
is an additional complication: a harmonic oscillator is 
a quantum system of infinite dimension, and no finite 
amount of measurement data will constitute a quorum. 
In order to achieve reconstruction, one needs to make 
certain assumptions that limit the number of free pa- 
rameters defining the state in question. For example, the 
filtered back-projection imposes low pass filtering onto 
the Fourier image of the Wigner function, i.e. assumes 
the ensemble to possess a certain amount of "classical- 
ity" (Vogel, 2000). Such smoothing reduces the accu- 
racy of the reconstruction (Herman, 1980; Leonhardt, 
1997) and introduces characteristic ripples (Brcitcnbach 
et al, 1997) on the reconstructed phase-space density 
[Fig. da)]. 



Although errors cannot be eliminated completely, we 
would like a reconstruction method that guarantees a 
physically plausible ensemble and minimizes artifacts. 
This requirement is satisfied by the Maximum Likelihood 
(MaxLik) approach, which aims to find, among the vari- 
ety of all possible density matrices, the one which maxi- 
mizes the probability of obtaining the given experimental 
data set and is physically plausible. Because this method 
is relatively new, but is rapidly gaining popularity, we 
here present its relatively detailed description. A yet 
more comprehensive review on quantum MaxLik (limited 
to the discrete domain) is given in Hradil et al. (2004). 

2. Classical algorithm 

We begin with a brief discussion of the classical 
expectation-maximization method. Consider a certain 
system characterized a set of parameters f (such that 
Ti > and X^j'^i ~ 1)' which we need to determine. 
We are allowed to subject the system to a measurement 
with a random outcome. The probability of each possible 
result (indexed by j) is related to r linearly: 

W?{j) ^^r.hij, (32) 

i 

where all hij are known positive numbers. The measure- 
ment is repeated N times, of which each outcome occurs 
fj times. The goal is to infer the parameter set r from 

the set of measurement results /. 

The ideal inference is the one that satisfies the system 
of linear equations 

^i^Y.'-^h^r (33) 

i 

However, a solution to this system exists only if the num- 
ber of parameters is larger than the number of equations. 
Otherwise, we have to settle for less: find the distribution 
r which would maximize the probability (likelihood) 

C{r) = Y[[Wr{j)]^' , (34) 

j 

of the observed measurement result. This approach has 
a very large variety of applications ranging form image 
de-blurring to investment portfolio optimization. 

The maximum-likelihood parameter set is determined 
by the so-called expectation-maximization (EM) algo- 
rithm which consists of sequential iterations (Dempster 
et al, 1977; Vardi and Lee, 1993): 

h 

initialized with some positive vector r. Each single iter- 
ation step is known to increase the likelihood. Further- 
more, because the likelihood is a convex function (i.e. 
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> 



for any two distributions ri and r2 holds £( 

'^^''^"'2^^^''^"' )' ^'^^ iterations will approach the global like- 
lihood maximum. 



3. The discrete quantum case 

A quantum tomographic procedure can be associated 
with a positive operator- valued measure (POVM), with 
each possible measurement result described by a positive 
operator T\j , which occurs with a probability 



prJj) = Tr[n,p]. 



(36) 



Here, again, we are dealing with a linear inversion prob- 
lem, because the probabilities are proportional to the 
density matrix elements. However, the latter are not nec- 
essarily positive (not even real) and their sum is not equal 
to one, so the EM algorithm in its original form has only 
limited application to the quantum case. 

In order to reconstruct a quantum state, we introduce 
the non-negative operator 



-Jj—f{. 



(37) 



As shown by Hradil (1997), the state that maximizes the 
likelihood ( 34 1 obeys the extremal equation 



as well as 



R{Pq)Pq = PoR{Po) = Pa, 



R{po)poR{po) = Po- 



rn 



(39) 



This iteration ensures positivity of the density matrix 
and has shown fast convergence in a variety of exper- 
iments. However, there is no guarantee for monotonic 
increase of the likelihood in every iteration; on the con- 
trary, there exists a (somewhat pathological) counterex- 
ample (Rehacek et ai, 2007). There remains a risk that 
the algorithm could fail for a particular experiment. 

The remedy against this risk is proposed in Hradil 
et al. (2004) and further elaborated in Rehacek et al. 
(2007). These publications present a "diluted" linear it- 
eration 



P 



(fc+i) 



1 + e ^ 1 + e 



(41) 



which depends on a single parameter e that determines 
the "length" of the step in the parameter space associated 
with one iteration. For e — > cxd, the iteration becomes 
RpR. On the other hand, in the limit of e ^ 0, there is a 
proof that the likelihood will monotonically increase and 
the iterations will converge to the maximum-likelihood 
state. We thus obtain a backup algorithm for the case 
the likelihood fails to increase in the RpR iteration. In 
practice, however, this situation is not likely. 

In some tomography schemes, one or more possible 
measurement results may not be accessible and, conse- 
quently, G = J2j is not equal to the unity operator. 



Then the extremal map ( 39 1 should be replaced by 



G-'R{po)poR{po)G-' = po 



(42) 



to avoid biased results (Rehacek et ai, 2001; Hradil et 
al, 2006; Mogilevtsev et ai, 2007). This issue may be- 
come significant in homodyne tomography reconstruc- 
tion, which we discuss next. 



One can intuitively understand these equations as fol- 
lows: when p is the maximum-likelihood state, we have 
fj/N « pr^, so the operator R becomes J^j flj, which is 
normally unity. 

The analogy to the classical scheme would suggest 
an iterative procedure = R{p'^^'^)p'^'''^ based on 

Eq. (38 1. However, unfortunately, such iteration does 



not preserve positivity of the density matrix [unless it 
is guaranteed to be diagonal in some basis, in which 
case the iteration reduces to Eq. ( pS] ) — such as in Ba- 
naszek (1998a, b)]. A possible solution is to apply the 
expectation-maximization iteration to the diagonalized 
density matrix followed by re-diagonalization (Rehacek 
et al, 2001; Artilles et al, 2005). 

A more common approach to constructing the iterative 
algorithm relies on Eq. ([39]) (Hradil et al, 2004). We 
choose some initial denstity matrix as, e.g. 
and apply repetitive iterations 



P 



(fc+i) 



-N i?(p('=))p('=)i?(/5('=)) 



(40) 



where denotes normalization to a unitary trace. Here- 
after we refer this scheme as the "RpR algorithm" . 



4. Iterative scheme for homodyne tomography 

The applications of MaxLik to homodyne tomography 
have been pioneered by Banaszek (1998a, b), who re- 
constructed the photon-number distribution (the diag- 
onal density matrix elements which correspond to a 
phase-randomized optical ensemble) from a Monte-Carlo 
simulated data set by means of the classical EM algo- 
rithm. This idea was then extended to reconstructing the 
Wigner function point-by-point (Banaszek, 1999) by ap- 
plying phase-dependent shifts to the experimental data. 
In a subsequent publication, Banaszek et al (1999) dis- 
cussed direct MaxLik estimation of the density matrix, 
but presented no specific algorithm. More recently, the 
RpR iterative algorithm was adapted to OHT (Lvovsky, 
2004) and has since been widely used in experiments on 
homodyne reconstruction. We describe this adaptation 
below. 

For a given local oscillator phase 6, the probability to 
detect a particular quadrature value Qg is proportional 
to 



prJQe,d)^Tr[tliQg,9)p], 



(43) 
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where Il{Qe,9) is the projector onto this quadrature 
eigenstate, expressed in the Fock basis as 



{m\fl{Qe,0) \n) - {m\Qe,d) {Qe,e\n) 



(44) 



where the wavefunction {mlQe 



„i7n0 



by Eq. (261 



i>miQe) is given 



Because a homodyne measurement generates a num- 
ber from a continuous range, one cannot apply the itera- 
tive scheme (40 1 directly to the experimental data. One 



way to deal with this difficulty is to discretize the data 
by binning it up according to 6 and Qg and counting 
the number of events /g^.e belonging to each bin. In 
this way, a number of histograms, which represent the 
marginal distributions of the desired Wigner function, 
can be constructed. They can then be used to imple- 
ment the iterative reconstruction procedure. 

However, discretization of continuous experimental 
data will inevitably lead to a loss of precision^. To lower 
this loss, one needs to reduce the size of each bin and 
increase the number of bins. In the limiting case of in- 
finitely small bins, fq^^e takes on the values of either 
or 1, so the likelihood of a data set {{Qi, Oi)} is given by 



(45) 



and the iteration operator (37 1 becomes 



(46) 



where i — 1 ... TV enumerates individual measurements. 



The iterative scheme ( 40 1 can now be applied to find the 



density matrix which maximizes the likelihood ( 45 1 



In practice, the iteration algorithm is executed with 
the density matrix in the photon number representation. 
In order to limit the number of unknown parameters, we 
truncate the Hilbert space by excluding Fock terms above 
a certain threshold. This is equivalent to assuming that 
the signal field intensity is limited. In many experimen- 
tal situations, application of this assumption is better 
justified than the low-pass filtering used in the filtered 
back-projection algorithm. 

Fig. |3] compares the inverse linear transform and Max- 
Lik reconstruction methods in application to the experi- 
mental data from Lvovsky and Mlynek (2002). The data 
set consists of 14152 quadrature samples of an ensem- 
ble approximating a coherent superposition of the single- 
photon and vacuum states. We see that the MaxLik 
method eliminates unphysical features and artefacts that 
are present in the inverse Radon reconstruction. 



In fact, recent research shows the precision loss due to binning to 
be insignificant. On the other hand, binning greatly reduces the 
number of data and thus expedites the iterative reconstruction 
algorithm (D. Mogilevtsev, 2007, private communication) 



5. Error handling 



A homodyne detector of non-unitary efficiency 77 can be 
modeled by a perfect detector preceded by an absorber. 
In transmission through this absorber, photons can be 
lost, and the optical state undergoes a so-called gener- 
alized Bernoulli transformation (Leonhardt, 1997). If 77 
is known, the Bernoulli transformation can be incorpo- 
rated into the matrices of the POVM elements Il{Qg, 9) 
(Banaszek et ai, 1999; Lvovsky, 2004). These operators 
can then be used to construct the matrix R, so the itera- 
tive algorithm will automatically yield the density matrix 
corrected for detector inefficiencies. 

Theoretically, it is also possible to correct for the detec- 
tor inefficiencies by applying the inverted Bernoulli trans- 
formation after an efficiency-uncorrected density matrix 
has been reconstructed (Kiss et ai, 1995). However, this 
may give rise to unphysically large density matrix ele- 
ments associated with high photon numbers. Similar con- 
cerns about possible numerical instability arise when the 
detector inefficiency is being accounted for in the pattern- 
function reconstruction (Kiss et ai, 1995). With the inef- 
ficiency correction incorporated, as described above, into 
the MaxLik reconstruction procedure, this issue does not 
arise (Banaszek, 1998b). 

Another source of error in OHT MaxLik estimation 
can be the incomplete character of the homodyne mea- 
surements: the sum of the projection operators G — 
J2i^iQiT^i) is not equal identity operator (even in the 
truncated Fock space). Mogilevtsev et al. (2007) found 
that the deviation can be quite significant. This issue 
can be resolved by employing the iteration based on the 



biased extremal equation (42| instead of Eq. (p59|, such 
as in the experimental worE 



by Fernholz et al. (2008) 



Finally, we discuss statistical uncertainties of the re- 
constructed density matrix. In generic MaxLik algo- 
rithms, they are typically estimated as an inverse of the 
Fisher information matrix (Rao et al., 1945; Cramer, 
1946). This method can be generalized to the quantum 
case (Hradil et al., 2004; Usami et al., 2003). In applica- 
tion to OHT, calculating the Fisher information appears 
quite complicated due to a very large number of indepen- 
dent measurements involved. 

A sensible alternative is offered by a clumsy, yet sim- 
ple and robust technique of simulating the quadrature 
data that would be associated with the estimated den- 
sity matrix pml if it were the true state. One generates 
a large number of random sets of homodyne data accord- 
ing to Eq. ( 43 1 , then applies the MaxLik reconstruction 



scheme to each set and obtains a series of density matri- 
ces p'f^, each of which approximates the original matrix 
Pml- The average difference {\pml — p'k\)k evaluates the 
statistical uncertainty associated with the reconstructed 
density matrix. 
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C. Maximum-entropy reconstruction 

The maximum-entropy (MaxEnt) method is appHed in 
the situation opposite to that of the MaxLik approach: 
when the number of equations in system (33) (i.e. the 



number of available data) is smaller than the number 
of unknown parameters (Buzek, 2004). In this case, the 
solution is not unique, and MaxEnt looks for the least 
biased solution, i.e. the one that maximizes the von Neu- 
mann entropy S = — Tr(/51og p). 

In OHT, one usually collects a very large (10^-10^) 
number of data points, hence the situation where the 
MaxEnt method is applicable is uncommon. Neverthe- 
less, Buzek and Drobny (2000) have elaborated appli- 
cation of this method to homodyne tomography and 
performed reconstruction of various simulated data sets. 
They found the results to be significantly better than 
those obtained by inverse linear transform, particularly 
in situations of incomplete tomographic data (marginal 
distributions available for a small number of phases or 
measured on short intervals). 



IV. TECHNICAL ASPECTS 

A. Time-domain homodyne detection 

When homodyne detection was first introduced to 
quantum optical measurements in mid-1980s, it was used 
for evaluating field quadrature noise rather than full 
state tomography. Such measurements are convenient 
to perform in the frequency domain, observing a cer- 
tain spectral component (usually around 1-10 MHz where 
the technical noise is minimized) of the photocurrent 
difference signal using an electronic spectral analyzer. 
Frequency-domain detection was used, for example, to 
observe quadrature squeezing (Slusher et al, 1985; Wu 
et al, 1986). 

Quantum-information applications require measure- 
ment of optical modes that are localized in time. Homo- 
dyning has to be performed in the time domain: differ- 
ence photocurrent is observed in real time and integrated 
over the desired temporal mode to obtain a single value 
of a field quadrature. Repeated measurements produce 
a quantum probability distribution associated with this 
quadrature. 

In this section, we discuss the design of time-domain 
balanced detectors that operate with pulsed local os- 
cillators. The first such detector was implemented by 
Smithey et al. (1992,1993a) in their original quantum 
tomography experiments. Among subsequent schemes 
we note that of Hansen et al. (2001) which features a 
higher bandwidth and a signal-to-noise ratio as well as 
that of Zavatta et al. (2002), exhibiting a further signifi- 
cant bandwidth increase at a cost of a somewhat poorer 
noise characteristics. 

Figure |4|^a) shows the main elements of the circuit of 
Hansen et al. (2001), which are typical for today's pulsed. 



time-domain HDs. A pair of high-efficiency photodiodes 
are wired in series to subtract their output currents, and 
this difference signal is amplified by a charge-sensitive 
transimpedance amplifier, followed by a pulse-forming 
network. The optics in front of the photodiodes per- 
mits variable attenuation of the input to each photodiode 
(alternatively, a setting with two polarizing beam split- 
ters and a half-wave plate between them can be used for 
combining the local oscillator and the signal). Thorough 
balancing of the photodiodes' photocurrents is essential 
for the proper operation of the BHD. 

With each LO pulse, the detector produces a burst of 
amplified subtraction photocurrent [Fig. |4|b)]. Because 
the response time of the detector is much slower than 
the width of the laser pulse, the generated signal is pro- 
portional to the time integral of the photocurrent over 
the pulse duration. This is a single sample of the field 
quadrature noise in the spatiotemporal optical mode of 
the local oscillator pulse. 

To prove that the pulsed noise generated by the homo- 
dyne detector with a vacuum signal input is indeed the 
shot noise, one needs to verify that the output rms noise 
scales as the square root of the LO power^ [Fig. |4|c)]. 
This is a signature distinguishing the shot noise from the 
classical noise (proportional to the local oscillator inten- 
sity) and the electronic noise (which is constant) (Bachor 
and Ralph, 2004). 

Design of time-domain BHD is more technically chal- 
lenging than its frequency-domain counterpart. First, 
the electronics must ensure time separation of responses 
to individual laser pulses. The shot-noise difference 
charge must be low-noise amplified within a bandwidth 
exceeding the local oscillator pulse repetition rate. Sec- 
ond, precise subtraction of photocurrent is necessary in 
order to eliminate the classical noise of the local oscil- 
lator. There is a competition between this requirement, 
which is easier satisfied at lower LO energies, and that 
of a sufficiently strong subtraction signal TV- , which in- 
creases with the LO power. The compromise is achieved 
on the scale of ~ lO'^-lO^ photoelectons. Finally, 
the measured quadrature values must not be inffuenced 
by low- frequency noises. The detector must thus provide 
ultra-low noise, high subtraction and ffat amplification 
profile in the entire frequency range from (almost) DC to 
at least the LO pulse repetition rate. 

A typical dilemma faced by a BHD designer is a trade- 
off between the signal-to-noise (more precisely, shot-to- 
electronic noise) ratio and the bandwidth (Raymer and 
Beck, 2004). An amplifier with a higher bandwidth 
usually exhibits poorer noise characteristics (Nicholson, 
1974; Radeka, 1988). An additional bandwidth limita- 
tion arises from the intrinsic capacitance of photodiodes, 
which may cause instability in the amplification circuit. 



^ This follows from Eq. | |l3| l, because la^l = \/Nlo a^nd Qo varies 
on the scale of 1. 
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FIG. 4: (a) Electro-optical scheme of the homodyne detector, 
(b) A superposition of multiple oscilloscope traces of the de- 
tector output. Each pulse produces a time-resolved quantum 
noise sample, (c) RMS peak amplitude of the noise pulses 
as a function of the LO power showing the expected square 
root power dependence up to the LO intensities of 3 x 10* 
photons per local oscillator pulse. Filled squares show the 
measured noise variances, open squares have the electronic 
noise background corresponding to 730 electrons/pulse sub- 
tracted. From Lvovsky et al. (2001). 



Technologically, this capacitance is determined by the 
thickness of the photodiode PIN junction; reducing this 
thickness compromises the quantum efficiency. A homo- 
dyne detector with a time resolution capable to accom- 
modate a typical repetition rate of a mode-locked, pulsed 
Ti:Sapphire laser (around 80 MHz) has been demon- 
strated by Zavatta et al. (2002, 2005b). 

Suppression of the homodyne detector electronic noise 
is important for quantum state reconstruction. As shown 
by Appel et al. (2007), presence of the noise leads to an 
equivalent optical loss of l/S, where S is the detector's 
signal-to-noise ratio. 

We briefly note that time-domain homodyne detec- 
tion finds its applications not only in quantum tomo- 
graphy, but also in other fields of quantum and classical 
technology. One example is shot-noise-limited absorp- 
tion measurements at subnanowatt power levels achiev- 
able thanks to the very low technical noise (Hood et al., 
2000). Another is ultrafast, ultrasensetive linear optical 
sampling for characterizing fiber optical systems (Dorrer 
et al., 2003). Time-domain homodyning is also an es- 
sential element of continuous-variable quantum cryptog- 
raphy (Silberhorn et al, 2002; Grosshans and Grangier, 
2002; Funk, 2004; Lodewyck et al, 2007). 



B. Matching the mode of the local oscillator 

1. The advanced wave 

In homodyne detection, the spatiotemporal optical 
mode to be measured is determined by that of the local 
oscillator. In this way, OHT provides indirect informa- 
tion on the modal structure of the signal field. This is 
useful for evaluating quantum optical information pro- 
cessing systems, which require that interacting optical 
qubits be prepared in identical, pure optical modes. On 
the other hand, achieving the mode matching between 
the local oscillator and the signal, or even preparing the 
signal state in a well-defined, pure spatiotemporal mode, 
can be challenging. In this section, we discuss the mode 
matching techniques, specializing to an important par- 
ticular case of the signal state being a heralded single 
photon. 

In order to prepare a heralded photon, a paramet- 
ric down-conversion (PDC) setup is pumped relatively 
weekly so it generates, on average, much less than a sin- 
gle photon pair per laser pulse (or the inverse PDC band- 
width). The two generated photons are separated into 
two emission channels according to their propagation di- 
rection, wavelength and/or polarization. Detection of a 
photon in one of the emission channels (labeled trigger 
or idler) causes the state of the photon pair to collapse, 
projecting the quantum state in the remaining (signal) 
channel into a single-photon state [Fig. [5|a)]. Proposed 
and tested experimentally in 1986 by Hong and Mandel 
(1986) as well as Grangier et al. (1986), this technique 
has become a workhorse for many quantum optics exper- 
iments. 

The biphoton is a complex entangled state with many 
parameters (spectrum, direction, polarization, etc.) of 
the two photons highly correlated: 



*(w5,tJ(,/Cs,fct) 1 



l^^ j:^^ duisdujtdksdkt, 
(47) 

where uj and k denote the frequencies and wavevectors 
of the signal and trigger photons. If the trigger photon 
is measured with any uncertainty in one of these param- 
eters, the signal photon will be prepared in a non-pure 
state 



Ps^TrtT{uJt,kt)\^s^){'^s 



(48) 



where T{LUt, kt) is the transmission function of the filters 
in the trigger channel defining the measurement uncer- 
tainty. 

Formation of the heralded mode is nicely illustrated 
by the heuristic concept of advanced waves proposed by 
Klyshko (1998a, b, c) and further advanced by Aichele et 
al (2002). According to this concept, the trigger photon 
detector is replaced with a fictitious light source, which, 
at the moment of detection, produces a classical incoher- 
ent electromagnetic wave traveling backwards in space 
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FIG. 5: Parametric down-conversion and the advanced wave 
model, (a) Preparation of single photons by conditional mea- 
surements on a biphoton state, (b) The Klyshko advanced 
wave model. The trigger detector is replaced with an inco- 
herent light source, which generates an incoherent advanced 
wave propagating backwards in space and time. Nonlinear 
interaction of this wave with the pump produces a difference- 
frequency pulse that mimics that of the conditionally pre- 
pared photon, (c) In an experiment, a laser beam, aligned for 
maximum transmission through all the filters, can model the 
advanced wave. From Aichele et al. (2002). 



and time [see Fig.|5];b)]. When propagating through the 
trigger channel filters, the advanced wave acquires some 
degree of spatiotemporal coherence, quantified by the 
filters' width. It then enters the down-conversion crys- 
tal and undergoes nonlinear interaction with the pump 
pulse, generating a difference-frequency pulse. This pulse 
turns out to be completely identical, in its modal char- 
acteristics, to the optical mode of the conditionally pre- 
pared single photon (Aichele et al., 2002), and is thus 
helpful for visualizing many of its properties. 

Suppose, for example, that the pump is pulsed (femto- 
or picosecond). Because the timing jitter of the photon 
counter event is typically on the order of a nanosecond, 
the exact moment of when the photon pair has been emit- 
ted is uncertain, so the advanced wave can be assumed 
continuous in time. The nonlinear interaction between 
the advanced wave and the pump is however restricted by 
the spatiotemporal window determined by the coherent 
pump pulse. If the latter is much narrower than the co- 
herence time and coherence width of the advanced wave, 
the difference-frequency pulse will be almost transform 
limited, both in the spatial and temporal dimensions. 

We conclude that narrow spatial and spectral filtering 
of the trigger photon can be used to obtain the signal 
photon in a pure spatiotemporal mode. To our knowl- 
edge, for the first time this matter has been investigated 
theoretically by Zukowski et al. (1995), independently 
by Rarity (1995), and later confirmed in a more detailed 
study by Ou (1997). Specifically in the context of OHT, 
theoretical treatment was given by Grosshans and Grang- 
ier (2001) as well as Aichele et al. (2002). 

Mode purity of the signal photon does not by itself 
guarantee its matching to the local oscillator. The ad- 
vanced wave model suggests the following experimental 



procedure for achieving this matching. Although the ad- 
vanced wave propagates backwards in space and time and 
is thus a purely imaginary object, it can be modeled by 
a forward-going alignment beam inserted into the trigger 
channel so that it overlaps spatially and temporally with 
the pump beam inside the crystal and passes through the 
optical filters [Fig.[5]^c)]. Nonlinear interaction of such an 
alignment beam with the pump wave will produce differ- 
ence frequency generation into a spatiotemporal mode 
similar (albeit no longer completely identical) to that of 
the conditionally prepared single photon. If one observes 
and optimizes the interference pattern between this wave 
and the local oscillator, one can be sure that, after block- 
ing the alignment beam, the mode of the signal photon 
will be matched to that of the local oscillator (Aichele et 
al, 2002). 

Controlling the spatial mode of the signal photon is 
simplified if a single-mode optical fiber is used as an 
optical filter instead of a pinhole arrangement. Such a 
filter automatically selects a pure spatial mode in the 
trigger channel, which transforms to a spatially pure sig- 
nal photon. It is also advantageous in terms of the pair 
production rate [Ourjoumtsev et al. (2006a), supporting 
material] Unfortunately, there is no similar arrangement 
possible for the spectral (temporal) mode matching. 



2. Decorrelatmg photons 

Reducing the spectral line width of the trigger filter 
will improve the mode purity of the heralded photons, 
but also reduce their production rate. This compromise 
would be avoided if we could arrange the PDC setup in 
such a way that the trigger and signal photons in the out- 
put of the down-converter are uncorrelated: the function 
^' in Eq. ( 47 1 could be written as 



(49) 



In this case, detection of any photon in the trigger chan- 
nel signifies that the signal photon has been emitted into 
a pure spatiotemporal mode defined by the function ipg. 

The first detailed theoretical inquiry into preparation 
of uncorrelated down-conversion spectra was made in 
Grice et al. (2001), based on general theoretical analy- 
sis of Keller and Rubin (1997). This theory was further 
elaborated in U'Ren et al. (2005, 2007). The configura- 
tion of the correlation function ^> depends primarily on 
the energy-conservation condition 



UJs+LL!t= LUp 

and the phase-matching condition 



hp. 



(50) 



(51) 



For any generated pair of photons with parameters 
{uJs,ks,ujtikt) there must exist a pump photon (up^kp) 
for which the above equations are satisfied. 
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Suppose the PDC occurs in an almost coUinear config- 
uration and the crystal is aligned so that Eqs. (50 1 and 



(51 1 simultaneously hold for the central pump frequency 
ujpQ and some signal and idler frequencies Wso ^^nd cotOi 
respectively. The frequency and the wavevector are con- 
nected through dispersion relations: dw — VgrA\k\, where 
Vgr is the wave's group velocity. Neglecting dispersion 
orders higher than one (analysis beyond this approxima- 
tion is made in U'Ren el ai, 2005), we cast Eq. (51 1 into 
the form 



Vqr.t 



(52) 



Equations (50 1 and (52 1, plotted in the {ujs,ujt) plane 
with Up = LOpo, form straight lines crossing at {uJsQ,ujto)- 
Considered more accurately, the lines defined by these 
equations are not infinitely narrow. This is because the 
pump is pulsed, so it contains photons not only at ujpo, 



but in a finite frequency range 



< 



Tr/Tp deter- 



mined by the inverse pump pulse width Tp. The phase- 
matching condition also has a tolerance: \ks + kt ^ kp\ ^ 
tt/Lc, where Lc is the crystal length that limits the region 
of nonlinear optical interaction. The down-conversion 
spectrum is thus determined by the overlap of two band- 
shaped areas in the {tUs^ut) plane [Fig. [6](a,b)]. 

As can be seen, the spectrum does not automatically 
uphold Eq. (49 1. However, by choosing the crystal length 



and other parameters of PDC, one can engineer the tilt 
angle and the width of the spectral region in which phase 
matching is satisfied, such that the overlap area can be 
expressed in the product form ( 49 1 , implying uncorre- 
lated photon spectra. The condition that has to be ful- 
filled takes the form (Grice et ai, 2001) 



1 



= -0.048L-' 



1 



1 



1 



1 

Vqr.t 



(53) 



where L is the crystal length and a is the pump spectrum 
width. 

There exist several theoretical proposals on shaping 
the phase matching region. Grice et al. (2001) calculate 
that Eq. (53 1 satisfies in a BBO crystal for degenerate 



coUinear type-II down-conversion if the pump wavelength 
is set to 800 nm [Fig.[6](c,d)]. They calculate a number of 
alternative down-conversion configurations with decorre- 
lated spectra. U'Ren et al. (2003) propose to implement 
PDC in a slightly non-coUinear configuration and impose 
additional restrictions onto the signal-idler spectrum by 
collecting only the photons emitted at certain angles. In 
theoretical works by Walton et al. (2003, 2004), down 
conversion occurs in a nonlinear waveguide, pumped al- 
most orthogonally to the guided direction. This also 
allows restricting the transverse components of the sig- 
nal and idler photon momenta, but without compromis- 
ing the pair production rate and with the possibility to 
choose the central wavelength of each photon. Torres et 
al. (2005) proposes to engineer the down-conversion spec- 
trum by employing a chirped, tilted pump and utilizing 
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FIG. 6: Shape of the biphoton correlation spectrum deter- 
mined by (a): the energy conservation condition (501, (b): 



generic phase matching ( 52 1 in the case of coUinear type II 



PDC in a BBO crystal pumped at 400 nm, (c): engineered 
phase-matching (coUinear type II PDC in a BBO crystal 
pumped at 800 nm), (d): cumulative effect of conditions (a) 
and (c). The insert in (d) shows the Schmidt decomposition of 
the biphoton spectrum "^(ujsjijJi) ~ V-^mUm(ti's)i'm(wi); 
with a proper combination of geometrical paremeters of the 
experiment the Schmidt decomposition contains only one 
term, i.e. the spectrum (d) is uncorrelated. Reproduced with 
permission from Grice et ai, (2001). 



the Poynting vector walk-off effect. Raymer et al. (2005) 
put forward the idea of placing the down-conversion crys- 
tal into a microcavity, whose line width is much narrower 
than that allowed by the energy conservation and phase 
matching conditions. This leads a spectrally uncorrelated 
biphoton in a spatial mode defined by the cavity. U'Ren 
et al. (2006, 2007) theoretically show that the group de- 
lays can be controlled by means of a periodic assembly 
(superlattice) of nonlinear crystals and birefringent spac- 
ers, and perform a proof-of-principle experiment to this 
effect, albeit without actually achieving an uncorrelated 
spectrum. 

The only experimental demonstration of a virtually un- 
correlated down-conversion spectrum to date is offered 
by Mosley et al. (2007). They use a relatively long 
potassium-dihydrogen-phosphate (KDP) crystal and a 
pump wavelength of 415 nm. Under these conditions, 
the pump will propagate with the same group velocity as 
the idler photon. Then the region allowed by the phase- 
matching condition [Fig. [6]jb)] becomes vertical and nar- 
row, so the overlap region exhibits almost no correlation. 
This is confirmed by observing high-visibility Hong-Ou- 
Mandel interference of heralded signal photons from sep- 
arate crystals. 
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FIG. 7: Spectrum of the Fabry-Perot parametric cavity. If 
one of the spectral modes (ut) is selected in the trigger chan- 
nel, a click of the trigger detector heralds production of the 
signal photon in a pure cavity mode at Us = Swo — i^t ■ 



3. The continuous-wave case 

A completely different approach to producing heralded 
photons must be taken if the pump is monochromatic 
(continuous) and down-conversion occurs in an optical 
cavity, such as in Neergaard-Nielscn et al. (2007). Spatial 
mode matching is simplified in this configuration because 
both photons are prepared in the spatial mode of the cav- 
ity. The biphoton spectrum is determined by the cavity 
transmission spectrum: it consists of narrow (a few MHz) 
equidistant modes separated by the cavity free spectral 
range (FSR)^. If the down-converter is pumped at a fre- 
quency 2luq (where luq coincides with one of the cavity 
resonances), down-converted photons will be generated 
at frequencies LUo±n x FSR (Fig. [7|. 

Under these circumstances, the time uncertainty of the 
trigger detector event is insignificant, so the advanced 
wave can be considered to be a short pulse. This elim- 
inates the need for narrow filtering of the trigger pho- 
ton; it is sufficient to apply a spectral filter that would 
transmit one of the cavity modes (e.g. ujt — ~ FSR). 
The advanced wave pulse will be filtered by the cavity. 
Its interaction with the pump will produce a difference- 
frequency pulse at ujs = ujQ + FSR which, in turn, gets 
filtered by the cavity. Because the cavity spectrum is ap- 
proximately Lorentzian, the temporal shape of the con- 
ditionally prepared mode is given by 



fit,tc) 



tcl 



(54) 



where tj, is the moment of the trigger event and 7 is the 
HWHM cavity hnewidth [Fig. [8];a)]9. 

Homodyne detection of the field in this mode requires 
a local oscillator pulse with the temporal shape [gi, {t) in 
Eq. ([s])] identical to f{t,tc). Such a pulse can be "tai- 
lored" from a continuous laser field by means of acousto- 



The free spectral range of a cavity equals the inverse roundtrip 
time of a photon inside the cavity. 
® Eq. j54| has to be modified if the weak pumping limit is not 
applicable (MoLmer, 2006; Sasaki and Suzuki, 2006). 



FIG. 8: Temporal modes of the trigger (1) and signal (2) 
photons. Assuming that the time resolution of the trigger 
detector is very high, the signal mode is defined by Eq. ( 54 1 . 



Cases (a) and (b) correspond to the weak and strong pump- 
ing regimes, respectively. Reproduced with permission from 
M0lmer (2006) 



or electro-optical amplitude modulation. An alternative, 
more practical, procedure, consists on using a continuous 
local oscillator {gL{i) = 1) and continuously acquiring 
the difference photocurrent as a function of time. The 
acquired photocurrent is then post-processed by multi- 
plying it by /(t, tc) and subsequently integrating over 
time. As evidenced by Eq. (10 1, the integrated differ- 



ence charge is the same as that obtained with a pulsed 
LO. This idea was first utilized in Neersgaard-Nielsen et 
al. (2006) and subsequently in Wakui et al. (2006) and 
Neersgaard-Nielsen et al. (2007). 



^. Strong pumpmg mode 

We now briefly discuss the regime of strong-pump 
PDC, such that the number of pairs generated within the 
time period corresponding to the inverse down-conversion 
bandwidth is not negligible. This case is complicated and 
largely uninvestigated for the following reason: if there 
exists any correlation in the joint time- or spectral distri- 
butions of the signal and trigger photons, narrow filtering 
of the trigger channel does not guarantee purity of the 
signal state. For example, in the monochromatic-pump 
case, a trigger event at time tc heralds the presence of 
a photon in the mode f(t,tc) but does not ensure that 
no incoherent contributions are present from photons in 
"nearby" modes f{t,t'c) [Fig. [sj^b)]. In order to prepare a 
high-purity state, one needs to have a high-efficiency trig- 
ger detector, which will not only trigger a pair production 
event, but ensure there are no more events nearby. If the 
state to be prepared is more complex than a Fock state, 
such as a the Schrodinger kitten state, the situation is 
even more complicated. 

The strong pumping case in the context of preparing 
single and multiple heralded photons was investigated by 
M0lmer (2006) as well as Nielsen and M0lmer (2007a,b). 
Sasaki and Suzuki (2006) report a comprehensive theo- 
retical study employing mode expansion in the basis of 
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prolate spherical functions and obtain analytical expres- 
sions for a few limiting cases. 

Experimental handling of the mode mismatch in 
the strong pumping case is, on the other hand, rela- 
tively straightforward (Ourjoumtsev et al., 2006a, 2006b, 
2007). One introduces an empiric probability ^ that the 
state heralded by a click in the trigger detector belongs 
to the mode analyzed by the homodyne detector. With 
probability 1 — the heralded state belongs to an orthog- 
onal mode, which is equivalent to a dark count event. 
The actual value of ^ can be found by matching to the 
experimental data statistics. 



V. APPLICATIONS IN QUANTUM 
TECHNOLOGY 

Implementation of light for the purposes of QI tech- 
nology relies on our ability to synthesize, manipulate, 
and characterize various quantum states of the elec- 
tromagnetic field. OHT is used in optical quantum 
information as a way to solve the last of the above 
tasks. In this section, we discuss applications of OHT 
to "discrete-variable" quantum-optical information and 
quantum-optical technology in general. We review the 
new states of light that have been created in the last 
few years, methods of their preparation, and their tomo- 
graphic reconstruction. 

It is convenient to restrict this review by the temporal 
boundaries of the present century. The only nonclassical 
state of light investigated by OHT prior to 2001 was the 
squeezed state (Smithey et al., 1993; Breitenbach et al., 
1997). The last few years, on the contrary, have shown a 
technology boom, resulting in a plethora of new quantum 
optical states (Table 1) , some of which are significant not 
only to QI technology, but to the very roots of quantum 
physics. 



A. Fock state tomography 

The first non-Gaussian state to be studied by OHT 
(Lvovsky et al., 2001) is the single photon. This is not 
surprising, given the role this state plays in basic and 
applied quantum optics. Another important motivation 
for this experiment was to demonstrate reconstruction of 
an optical state whose Wigner function takes on negative 
values. 

The schematic and results of the experiment are shown 
in Fig. [9] The experiment employed a picosecond 
Ti:sapphire laser at a 790-nm wavelength. Pulsed single 
photons were prepared by conditional measurements on a 
biphoton state generated via parametric down-conversion 
(in the weak pumping regime). Narrow spatiotemporal 
filtering of the trigger photon was used as outlined in 
Sec. |IVB[ The field state in the signal channel was char- 
acterized by means of optical homodyne tomography. 



Reference 

Lvovsky et al. (2001) 
Zavatta et al. (2004b) 
Neergaard-Nielsen (2007) 
Lvovsky and Babichev (2002) 
Lvovsky and Mlynek (2002) 
Babichev et al. (2003) 
Babichev et al. (2004a) 
Zavatta et al. (2006a) 
Babichev et al. (2004b) 
Zavatta et al. (2006a,b) 
Zavatta et al. (2004a) 
Zavatta et al. (2005a) 
Zavatta et al. (2007) 

Parigi et al. (2007) 

Ourjoumtsev et al. (2006b) 
Wenger et al. (2004b) 
Neergaard-Nielsen et al. (2006) 
Ourjoumtsev et al. (2006a) 
Wakui et al. (2007) 
Ourjoumtsev et al. (2007) 



State 

Single-photon Fock state 
|1) 



Displaced single photon 
Single-rail qubit 
q|0)+/3|1) 



Dual-rail qubit 

q|0,1)+/3|1,0) 

Photon-added coherent state 

at|a) 

Photon-added thermal state 



Two-photon Fock state |2) 
Photon-subtracted squeezed 
state ( "Schrodinger kitten") 
1) + a |3), where \a\ < 1 



Squeezed "Schrodinger cat" 



TABLE I: Quantum states recently characterized by OHT 

Remarkably, all imperfections of the experiment (losses 
in transmission of the signal photon, quantum efficiency 
of the HD, trigger dark counts, mode matching of the 
signal photon and the local oscillator, and spatiotemporal 
coherence of the signal photon) had a similar effect on the 
reconstructed state: admixture of the vacuum |0) to the 
ideal Fock state |1): 

Pmcas = ?7|l)(l| + (l-^/)|0)(0|. (55) 

The greater the efficiency rj, the deeper the "well" in the 
Wigner function; classically impossible negative values 
are obtained when rj > 0.5. The original 2001 experiment 
showed rj = 0.55 ±0.01; later, this value was improved to 
0.62. 

An interesting feature of the optical single-photon state 
reconstruction is that the technique of homodyne tomo- 
graphy can be fully understood in the framework of clas- 
sical physics. This measurement could have been con- 
ducted (and interpreted) by someone who does not be- 
lieve in quantum mechanics. Yet the result of negative 
quasi-probabilities would appear absurd, incompatible 
with classical physics — thus providing a very strong 
evidence of "quantumness" of our world. 

A version of the single-photon Fock state tomography 
experiment, featuring a fast homodyne detector, allowing 
measurements at a full repetition rate of the pulsed laser 
(82 MHz), was reported by Zavatta et al. (2004b). 

Tomography of the two-photon Fock state |2) was re- 
ported by Ourjoumtsev et al. (2006b). The experimental 
arrangement is similar to the of Fig.|9ja), but the trigger 
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FIG. 9: The experiment on quantum tomography of the 
single-photon Fock state, (a) Simplified scheme; (b) 45000 
raw quadrature noise samples for the vacuum state and the 
Fock state; (c) the density matrix (diagonal elements) re- 
constructed using the quantum state sampling method; (d) 
the reconstructed Wigner function is negative near the origin 
point because the measurement efficiency reaches 62%. Side 
projections show phase-randomized marginal distributions for 
the measured vacuum and Fock states. From Lvovsky et 
al. (2001); Lvovsky and Babichev (2002). 



channel is split and directed into two single-photon de- 
tectors, whose simultaneous click triggers a homodyne 
measurement. In order to obtain a sufficient rate of 
such events, the parametric gain has to be non-negligible, 
which significantly complicates the analysis of the exper- 
iment. Ourjoumtsev and co-workers found that the ex- 
perimentally observed state can be fit by a theory taking 
into account five experimental parameters: 

• gain of the down-converter; 

• excess gain of a fictitious phase-independent ampli- 
fier placed behind the (ideal) down-converter; 

• BHD efficiency (including optical losses); 

• BHD electronic noise (whose effect is identical to 
optical loss, see Sec. |IV A l; 



• probability ^ that the single-photon detection 
events corresponds to a heralded mode that 



matches the local oscillator (see Sec. IV B I 



By optimizing these parameters, negative values of the 
experimental Wigner function were obtained. 



B. The optical qubit 

1. The dual-rail qubit 



As discussed in Sec. IBl one application of OHT, 



where it can be of advantage compared to other state 
characterization methods, is measurement of systems of 
dual-rail optical qubits. Tomography of one dual-rail 
qubit was performed by Babichev et al. (2004b) A 
dual-rail qubit, described by the state 



I* 



dual- 



-rail) =t|U,Ob) -p|0^,1b) 



(56) 



is generated when a single photon |1) incident upon a 
beam splitter with transmission t'^ and reflectivity p^, 
entangles itself with the vacuum state |0) present in the 
other beam splitter input. To perform tomography mea- 
surements, BHDs (associated with fictitious observers Al- 
ice and Bob) were placed into each beam splitter out- 
put channel [Fig. 10 'a)]. With every incoming photon, 
both detectors made measurements of field quadratures 
QeA and QeA with the local oscillators' phases set to 
9a and 9b, respectively. Fig. [To|b) shows histograms 
of these measurements, which are the marginal distribu- 
tions of the four-dimensional Wigner function of the dual- 
rail state. They have been used to determine the state via 
the maximum- likelihood technique (Sec. IIIB), resulting 

As expected. 



10 'c 



in the density matrix shown in Fig. 
the reconstruction reveals all the terms in the density 
matrix, including those (e.g. double-vacuum |0, 0)(0,0|) 
usually missed by the photon-counting method. 

Time-encoded dual-rail optical qubits were prepared 
and characterized by Zavatta et al. (2006a). The trig- 
ger channel of the PDC entered a fiber Michaelson in- 
terferometer whose path length difference was equal to 
the optical path inside the cavity of the master mode- 
locked Ti:Sapphire laser. The trigger detector, placed at 
the output of the interferometer, is then unable to dis- 
tinguish between a photon generated by some nth pump 
pulse that has traveled the long path of the interferom- 
eter and a photon generated by the {n -\- l)th pulse that 
has taken the shorter path. The conditional state was 
thus prepared in a coherent superposition 



^(n) Q(n+1) 



o(»)^i( 



n+l 



(57) 



where the phase (j) could be controlled by one of the in- 
terferometer mirrors. For characterizing the above state, 
a single homodyne detector suffices, but its signal has to 
be acquired at two different moments in time. 



A deterministic OHT scheme for two-mode state reconstruction 
in the Fock basis was first proposed and tested numerically in 
Raymer et al. (1996). A detailed theoretical analysis of different 
aspects of such an experiment was also made by Jacobs and 
Knight (1996) as well as Grice and Walmsley (1996). 
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a) jSi^JCi b) 86 = 59 = 71/2 86 = 71 
LOs 




FIG. 10: The experiment on homodyne tomography of 
the dual-rail qubit. (a) Scheme of the experimental setup, 
(b) Histograms of the experimental quadrature statistics 
pr^g(Q9A, Qes) for a symmetric beam splitter. Phase- 
dependent quadrature correlations are a consequence of the 
entangled nature of the state l^qubit). Also shown are indi- 
vidual histograms of the data measured by Alice and Bob, 
which are phase-independent, (c) Density matrix (absolute 
values) of the measured ensemble in the photon number rep- 
resentation. From Babichev et al. (2004b) 



2. Nonlocality of the single photon 



Whether the state (56) can be considered entangled 
is a widely debated issue. This controversy seems to be 
related to the wave-particle duality of light. If the pho- 
ton is viewed as a state of the electromagnetic oscillator, 
the notation ( |56[ ) is valid and denotes an entangled en- 
tity (Van Enk, 2005). If, on the other hand a photon is 
considered to be a particle, i.e. not a state but a carrier 
of a state, e.g. of a polarization state, the dual rate qubit 
should be written as a superposition of two localizations 
of one photon, which may not be seen as entangled. Ad- 
vocates of the former view proposed experiments on using 
the split single photon to demonstrate quantum nonlocal- 
ity (Oliver and Stroud, 1989; Tan et al., 1991; Banaszek 
and Wodkiewicz, 1999; Jacobs and Knight, 1996; Hessmo 
et al, 2004); others disputed them (Greenberger et al, 
1995; Vaidman, 1995). 

It is in the inherent nature of OHT to interpret the 
photon as a state of a field rather than a particle "in its 
own right" . In experiments on homodyne tomography 
of the delocalized photon, Babichev et al. (2004b) and 
Zavatta et al. (2006b) present different arguments that 
OHT characterization of the dual-rail qubit ( 56 ) can be 



means of a fictitious discriminator. Correlations between 
the discriminator outputs acquired by Alice and Bob ex- 
hibited a Bell- like interference pattern. For sufficiently 
high threshold values, its amplitude exceeds 1/a/2 and 
the Bell inequality is violated. Zavatta et al. (2006b) de- 
termined the Wigner function of the measured dual-rail 
state and showed it to violate the Bell test of Banaszek 
and Wodkiewicz (1999). 

Further evidence of the entangled nature of the de- 
localized single photon is its applicability as a resource 
in quantum communication protocols such as quantum 
teleportation and remote state preparation (RSP)^^. 



3. Remote state preparation using the nonlocal single 
photon state 

To implement RSP, Alice performs a measurement on 
her share of the entangled resource in a basis chosen in 
accordance with the state she wishes to prepare. De- 
pendent on the result of her measurement, the entangled 
ensemble collapses either onto the desired state at the 
receiver (Bob's) location, or can be converted into it by 
a local unitary operation. 

The experiment on tomography of the dual-rail qubit 
can be interpreted as an implementation of the remote 
preparation protocol in the continuous basis (Babichev et 
al, 2004a, Zavatta et al., 2006a). By performing a ho- 
modyne measurement on her part of the entangled state 
( 56 1 and detecting a particular quadrature value Qqa at 



the local oscillator phase 9, Alice projects the entangled 
resource (56 1 onto a quadrature eigenstate {QeA\- 



\M = {QeA^dAl"^) (58) 
= T {QeA, eA\l)A \0)b - P {QoA, Oa\0)a |1>b , 

which is just a coherent superposition of the single- 
photon and vacuum states, i.e. a single-rail optical qubit. 
By choosing her LO phase 9a and postselecting a partic- 
ular value of QeA, Alice can control the coefRcients in the 
superposition, i.e. remotely prepare any arbitrary state 
within the single-rail qubit subspace. 



4-. Teleportation using the nonlocal single photon state 

Entanglement contained in the delocalized single pho- 
ton state is between the single-photon and vacuum states. 
It allows rudimentary teleportation of single-rail qubits 



interpreted to violate Bell-type inequalities, albeit with 
loopholes. In Babichev et al. (2004), quadrature mea- 
surements have been converted to a dichotomic format by 



11 Both teleportation (Bennett et al, 1993) and RSP (Lo, 2000) are 
quantum communication protocols allowing disembodied trans- 
fer of quantum information between two distant parties by means 
of a shared entangled resource and a classical channel. The dif- 
ference between them is that in teleportation, the sender (Alice) 
possesses one copy of the source state, while in RSP she is instead 
aware of its full classical description. 
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|'0singie-raii) = ct |0)+/3 |1) by Hieans of a modified Bennet 
et al. (1993) protocol. Alice performs a Bell-state mea- 
surement on the source state and her share of |\E'duai-raii) 
by overlapping them on a beam splitter and sending both 
beam splitter outputs to single-photon detectors. If one 
of these detectors registers the vacuum state, and the 
other detects one photon, the input state of the Bell-state 
analyzer is projected onto |5'duai-raii) and Bob's channel 
obtains a copy of the source state (Pegg et al, 1998; 
Ozdemir et al., 2002). If the input state contains terms 
outside of the single-rail qubit subspace, these terms will 
be removed from the teleported ensemble. This is known 
as the "quantum scissors" effect. 

Although the implementation of the protocol requires 
highly-efficient, number- resolving photon detectors, its 
conceptual demonstration can be done with standard 
commercial units. Babichev et al. (2003) performed this 
experiment by using a weak pulsed coherent state as the 
source. The teleported state was characterized by means 
of OHT. The teleportation fidelity approached unity for 
low input state amplitudes but with this parameter in- 
creasing, it quickly fell off due to the effect of "quantum 
scissors" . 



5. Quantum- optical catalysis 




FIG. 11: The "quantum-optical catalysis" experiment, (a) 
The scheme. Measurements by the HD are conditioned on the 
single-photon detector registering a photon; (b) 14153 raw 
quadrature data; (c) absolute values of the density matrix 
elements in the Fock representation for a ~ 0.3. The beam 
splitter transmission is t^ — 0.075. From Lvovsky and Mlynek 
(2002). 



The two subsections above demonstrated how a single- 
rail optical qubit can be prepared by conditional mea- 
surements and linear-optical operations on a single pho- 
tons. Another way of achieving the same goal was re- 
ported by Lvovsky and Mlynek (2002). A single-photon 
state |1) and a coherent state \a) were overlapped on a 
high-reflection beam splitter. One of the beam splitter 
outputs was subjected to a measurement via a single- 
photon detector [Fig.[TT|a)]. In the event of a "chck" , the 
other beam splitter output is projected onto a single-rail 
qubit t|0) -t-a|l), being the beam splitter transmission 
[Figs. |3] and [ll];b,c)]. 



This result is somewhat counter-intuitive: a classi- 
cal coherent input state is converted into a nonclassi- 
cal single-rail qubit even though the input single pho- 
ton emerges "intact" at the output (hence the name 
quantum- optical catalysis). This transformation is an ex- 
ample of optical nonlinearity induced by a conditional 
optical measurement, the key principle behind linear- 
optical quantum computation (Knill et al., 2001; Koashi 
et al, 2001; Kok et al, 2007). 

We note that if the signal channel is analyzed without 
conditioning on the single-photon detection event, it ap- 
proximates another important nonclassical state of light, 
the displaced Fock state (Lvovsky and Babichev, 2002, 
and references therein). 



C. "Schrodinger cats" and "kittens" 



The Schrodinger cat (Schrodinger, 1935) is a famous 
Gedankenexperiment in quantum physics, in which a 
macroscopic object is prepared in a coherent superposi- 
tion of two classically distinguishable states. It brightly 
illustrates one of the most fundamental questions of 
quantum mechanics: at which degree of complexity does 
a quantum superposition of two states stops being a 
superposition and probabilistically becomes one of its 
terms? 

In quantum optics, the Schrodinger cat usually means 
a coherent superposition |a) ± |— a) of coherent states of 
relatively large amplitude and opposite phase (Buzek and 
Knight, 1995). In addition to the above fundamental as- 
pect, these states are useful for many quantum informa- 
tion protocols such as quantum teleportation (van Enk 
and Hirota, 2001), quantum computation (Ralph, 2003), 
and error correction (Cochrane et al, 1999). It is thus 
not surprising that experimental synthesis of Schrodinger 
cats has been an object of aspiration for several gener- 
ations of physicists. Recent years have marked a break- 
through: invention and experimental realization of two 
schemes that permit preparation of optical Schrodinger 
cats of arbitrarily high amplitudes. 
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1. Squeezing of single photons 

The first scheme was proposed by Dakna et al. (1997). 
An odd Schrodinger cat state of low amphtude can be 
decomposed into the Fock basis as follows: 



\a)-\-a) (X a|l) + — |3) + 



(59) 



For a ^ 1, this state is approximated, with very high fi- 
delity (Lund et al., 2004), by the squeezed single-photon 
state. Experimentally this state can be obtained by re- 
moving one photon from the squeezed vacuum 



|0) 



1 



V2 



C|2) 



e\4) 



(60) 



where C is the small squeezing parameter (Dakna et al, 
1997). The photon removal procedure consists of trans- 
mitting the state through a low-reflection beam splitter 
and sending the reflected mode to a single-photon detec- 
tor [Fig. 12 'a)]. A click in this detector indicates that at 
least one photon has been removed from |^'s). Because 
the reflectivity of the beam splitter is small, the prob- 
ability of removing more than one photon is negligible. 
The conditional state in the transmitted channel of the 
beam splitter is then approximated by 



*cond) OC -^C|l) 



Vsc^ 1 3) 



(61) 



Setting a — 6^ yields the "Schrodinger kitten" ( 59 1 



The first experiment implementing of this protocol was 
performed by Wenger et al. (2004b) and later improved 
by Ourjoumtsev et al. (2006a). Parametric deamplifi- 
cation of 150-fs, 40-nJ pulses at 850 nm in a 100 /im 
non-critically phase matched KNbOs crystal was used to 
generate pulsed squeezed vacuum (Wenger et al., 2004a). 
The heralded beam splitter output state was subjected to 
pulsed OHT, showing preparation of a Schrodinger kit- 
ten of size — 0.79 with a 70-% fidelity. Experimental 
imperfections could be modeled by the same five param- 
eters as in Ourjoumtsev et al. (2006b) (see Sec. VAl. 



In the continuous- wave regime, a similar procedure has 
been independently demonstrated by Neergaard-Nielsen 
et al. (20 06) [the results of this experiment are shown in 
Fig. ^h)] and by Wakui et al. (2007). 

These experiments, which combine, for the first time, 
the techniques of conditional preparation of single pho- 
tons and pulsed squeezing, are a significant technologi- 
cal breakthrough in quantum optical information tech- 
nology. However, the photon-subtracted squeezed state 
resembles a cat state only for small a's. An interesting 
method for generating Schrodinger cats of larger ampli- 
tudes has been proposed by Lund et al. (2004). Two 
small odd cat states"'^^ |CSS(a)), which we assume to be 



Under "even" and "odd" cat states we understand superpositions 
\q) + and \a) — \—q), respectively. 
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FIG. 12: (a) A generic scheme ol conditional preparation of 
the squeezed single-photon state ("Schrodinger kitten"), (b) 
The Wigner function and the density matrix reconstructed 
from the experimental data. Part (b) is reproduced with per- 
mission from Neergaard-Nielsen et al. (2006). 
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FIG. 13: A schematic for amplification of the "Schrodinger 
cat" state. The amplified cat emerges in channel g if detectors 
A and B click in coincidence. In the text, a — [3 — 7/\/2 is as- 
sumed. Reproduced with permission from Lund et al. (2004) 
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of equal amplitudes (a = overlap on a symmetric 
beam splitter BSl (Fig. |13| , which transforms them into 
an entangled superposition 



|CSS(a))JCSS(«)), 
\/2a, 



\/2c 



\/2c 



\/2a 



(62) 



of a cat state of amplitude v^a in one of the output 
channels and vacuum in the other. Now if we measure 
the state in channel g and find it to be not the cat state, 
the channel / will be projected onto the cat state. Such 
conditional measurement is implemented by overlapping 
channel g with a coherent state of amplitude 7 = ^/2a on 
an additional beam splitter BS2. If channel g contains 
a cat state, the interference will cause all optical energy 
to emerge at only one side of BS2. Therefore, detect- 
ing coincident photons in both outputs of BS2 indicates 
that channel g contained vacuum, and thus channel / is 
prepared in the Schrodinger cat state of amplitude %/2a. 

By applying this linear optical protocol repeatedly, we 
can "breed" Schrodinger cat states of arbitrarily high am- 
plitude. A remarkable practical advantage of this tech- 
nique is that it requires neither null single photon detec- 
tion nor photon number discrimination. It does however 
require a high degree of mode matching among the inter- 
fering optical channels; otherwise the fidelity will rapidly 
decrease (Rohde and Lund, 2007). 



2. Generating "Schrodinger cats" from Fock states 

An alternative way of generating optical cat states was 
recently proposed and implemented by Ourjoumtsev et 
al. (2007). This technique employs Fock states, rather 
than the squeezed state, as the primary resource. The 
procedure is remarkably simple: an n-photon number 
state is split on a symmetric beam splitter, and one of 
the output channels is subjected to homodyne detection. 
Conditioned on this measurement producing approxi- 
mately zero, the other beam splitter output mode will 
contain an approximation of the squeezed Schrodinger 
/n [Fig. 



14 'a 



cat of amplitude a 

To gain some insight into this method, let us assume 
that the quadrature measured by the preparation homo- 
dyne detector is the momentum P. The wavefunction of 
the initial Fock state is given by 



V'„(F) = {n\P) = H,,iP)exp 



p2 

' 2 



(63) 



In writing the above equation, we used Eq. (26), ne 



glected normalization factors and remembered that the 
Fock state is phase-independent (i.e. its wave function 
is the same for all quadratures). "Splitting" the state 
\n) means entangling it with the vacuum, which has the 
wavefunction Tpo{Po) = exp(— Pq/2), via transformation 
P ^ {P-Po)/V2, Po ^ {P + Po)/V2. Accordingly, 



the two-mode wavefunction of the beam splitter output 
is given by 



= i?«(^^)exp(^ — 



(64) 



Detecting the momentum quadrature value Po = in 
one of the modes has an effect similar to that discussed 
in Sec. |VB 3| it "remotely" prepares the other mode 
in the state with wavefunction 4>cond{P) ~ Hn{{P — 
Po)/V2)e-P'/\ 

This state is easier to analyze in the position quadra- 
ture representation. Making a Fourier transform of 

second 

(P), we find 



<^co„d(Q) = Q"e 



(65) 



This function has two peaks at Q — and vanishes 

at Q = and Q = ±00. The wave function of the coher- 
ent state |a), on the other hand, has a single maximum 
at q;-\/2- Thus the wavefunction (/)cat(a)(Q) of the even 
Schrodinger cat momentum-squeezed by a factor S has 
two peaks located at Q = ±SaV2. Good matching be- 
tween 0cond(Q) and 0cat(Q)(Q) obtains when their peaks 
have the same position and the same width, which hap- 
pens when a = y/n, S = 1/^/2- Amazingly, the fidelity 
of this matching increases with n, reaching the value of 
99 % already a,t n — 3. 

Experimental implementation of this protocol with 
n = 2 is almost identical to tomography of the two- 
photon Fock state (Ourjoumtsev et al., 2006b), except 
that two homodyne detectors are required. Reconstruc- 
tion of the output state bears close resemblance to the 
squeezed Schrodinger cat [Fig. [l4|b)], with differences 
mainly caused by technical limitations, which are similar 
to those in previous experiments by this group. 



D. Photon-added states 

Photon-added states (Agarwal and Tara, 1991) are gen- 
erated when the photon creation operator acts on an arbi- 
trary state [■(/') of light: \tp,Tn) = (a^)™ \ip). These states 
are nonclassical due to a vanishing probability of finding 
n < m photons (Lee, 1995). Recently, two important 
photon-added states were experimentally generated and 
characterized via OHT. 



1. Single-photon- added coherent states 

These states are of interest because in the limit of large 
a, they approximate highly classical coherent states |a) 
while for a — s- they become highly nonclassical Fock 
states |m). Therefore, photon-added coherent states can 
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a) 



FIG. 14: (a) A scheme of the experiment on conditional prepa- 
ration of the "Schrodinger cat" state from photon number 
states, (b) The experimental Wigner function. Reproduced 
with permission from Ourjoumtsev et al. (2007). 
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FIG. 15: (a) Conditional preparation of the single-photon- 
added coherent state, (b) With increasing a, the Wigner 
function of the reconstructed SPACS gradually evolves from 
a highly nonclassical to a highly classical shape. For \a\ = 2.6, 
also the Wigner function of the unexcited seed coherent state 
is shown. Reproduced with permission from Zavatta et al., 
(2004a). 



input coherent state results in gradual transition from the 
Fock state to an approximation of a coherent state. An 
interesting feature observed in SPACS of moderate am- 
plitudes is quadrature squeezing (up to 15%) associated 
with certain phases. 



be interpreted as a link between the particle and wave 
aspect of the electromagnetic field. 

Experimentally, photon addition can be implemented 
using a procedure opposite to photon subtraction de- 
scribed in the previous section. Instead of passing 
through a beam splitter, the target state is transmitted 
through a signal channel of a parametric down-conversion 
setup [Fig. [Tsjja)]. If a photon pair is generated in the 
down-converter, a photon is added to the target state. 
This event, which is heralded by a single photon emerg- 
ing in the trigger channel, can be followed by an OHT 
measurement of the signal ensemble. 

This scheme was implemented, for the first time, by 
Zavatta et al. (2004a, 2005a, 2005b). Thanks to high- 
bandwidth time-domain homodyne detection (Zavatta 
et al., 2002), no pulse-picking was necessary so the 
setup could be made highly compact and phase stable. 
Fig. |15[b) demonstrates how increasing amplitude of the 



2. Single-photon- added thermal states 

The thermal state is a phase-symmetric ensemble with 
Bose-Einstein photon number statistics. By itself, it is a 
classical state, but shows a high degree of nonclassicality 
when acted upon by the photon creation operator. This 
was shown experimentally by Zavatta et al. (2007). The 
thermal state was simulated by sending a coherent laser 
beam through a rotating ground glass disk and collecting 
a fraction of the scattered light with a single-mode fiber. 
It was then subjected to photon addition as described 
above, and subsequently to homodyne tomography. The 
measured state is verified to be highly nonclassical ac- 
cording to several criteria. 

In a related work, Parigi et al. (2007) apply a sequence 
of photon addition and subtraction operators to the ther- 
mal state. They find, contrary to classical intuition, but 
in full agreement with quantum physics, that the effects 
of these operators do not cancel each other, and, further- 
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more, depend on the sequence in which they are applied. 
This provides a direct evidence of noncommutativity of 
these operators, which is one of the cardinal concepts of 
quantum mechanics. 



VI. SPATIAL QUANTUM-STATE 
TOMOGRAPHY 



A. Spatial mode of the one-photon field 

As we discussed, homodyne tomography typically is 
used to measure the quantum state of light occupying a 
single selected optical mode, which is defined by the local 
oscillator pulse. If our goal is to characterize the field 
state in multiple modes, homodyne tomography becomes 
increasingly difficult. 

There is a special multi-mode situation, on the other 
hand, that is amenable to full characterization: If it is 
known a priori that only one photon (elementary exci- 
tation) of the field exists in a certain space-time volume, 
it is sensible to ask what is the temporal-spatial wave- 
packet mode that describes this photon. This task is 
close to that of finding the wave function of the photon 
treated as a massless particle. This notion is known to be 
controversial (see Smith and Raymer, 2007, for a review). 
However, if one restricts attention to the photon's trans- 
verse degrees of freedom in the paraxial approximation, 
the subtleties that arise can be circumvented. 

Assuming a constant polarization, a single-photon 
state of the quantized field can be represented by a su- 



perposition [c/. Eq. (Ill] 



I la) = a) I vac) 



d^kC{k)bl |vac) 



(66) 



= J d^kC{k) |1^), 



where 



\h)^\h)® n K^Cfc^') 



(67) 



is a one-photon state occupying a plane-wave mode with 
definite momentum p = fik. The function C{k) defines 
the spatial mode of the photon in the momentum repre- 
sentation. In the position representation for free space 
propagation, the matrix clement 



'■J+) 

E{r,t) = (vac|£; (f,t)|la) 
ie I £k 



(68) 



C{k) exp{ikf— iujrt), 



(2^)3£0 



defines the spatial distribution of the photon's field. The 
goal of spatial QST is to reconstruct the state by deter- 
mining the function E{f,t). 

Suppose a photon is created with a narrowly defined 
frequency wq and is propagating along the z axis with the 



wavenumber ko = loq/c. In the paraxial approximation, 
k^ kx , ky and 

fc.^y^r^-fco-|^-|^, (69) 
so we can rewrite Eq. (|68|) as 



E{r,t) = ieJ^^^^eM-i{^ot-koz)]u{r). (70) 
in which we define the spatial mode [with x = [x, y),kx = 

{kx : ky)] 



u{'r) ~ u{x, z) ^ / C{kx)exp 



2kn 



We compare the above expression with the Schrodinger 
evolution of a free particle of mass m in two dimensions, 
initially in a superposition ["(/'(O)) = J dk^ C{kx) \kx)'. 



il){x, t)^ C{kx) exp 



d?kx, (72) 



2m 



where m is the mass of the particle. Expressions (|71| and 
( [72| become equivalent if one sets z = ct and m = hko j c. 

We can utilize this equivalence by applying the pro- 
gram set forth in Eqs. ([T|--(|4]) in the Introduction in or- 
der to determine the transverse wave function of a pho- 
ton. We measure the beam intensity profile I{x, z) — 
(|i?(x, z)P) in different planes along the beam propa- 
gation direction. The transverse degrees of freedom of 
the wave evolve during propagation, allowing inversion of 
measured intensity (probability) distributions using the 
propagator ([s]), which now takes the form 



G{x\ x; t) oc exp 



ik{i\x 



2z 



(73) 



to determine the transverse wavefunction u{x, z) of the 
photon. 

If the transverse state of the photon is not pure, it is 
defined by the density matrix 



p{xi,X2) = {E{xi)E*{x2)) 



(74) 



with the angle brackets implying an ensemble average 
over all statistical realizations of the photon wave func- 
tion. Here we notice that the above definition is com- 
pletely analogous to that of the classical field correlation 
function determining the degree of its spatiotemporal co- 
herence^"^. Therefore, the tomography procedure we have 



13 In Sec. |IVB| we discussed the identity between the mode of 
the conditionally prepared photon and the classical difference- 
frequency signal generated by the advanced wave. This identity 
is expressed in terms of definition l |74| valid for both single pho- 
tons and classical fields. 
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developed for single photons is also applicable to classical 
fields, making them a useful "testing ground" for single- 
photon QST procedures. 

One can also introduce the transverse, two-dimensional 
spatial Wigner distribution at a particular plane in the 
fashion analogous to Eq.(16l: 



1 

4^ 



1 



1 



p{x+^Lx--£.)e~^''^<de (75) 
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where k,^. is the transverse-spatial wave-vector compo- 
nent. The transverse Wigner function is reconstructed 
from a set of beam intensity profiles using the in- 
verse Radon transform. Such a phase-space-tomography 
scheme was proposed in Raymer et al. (1994) for quan- 
tum or classical waves and implemented for the trans- 
verse spatial mode of a "classical" (coherent-state) light 
beam by McAHster et al. (1995). 



B. Non-interferometric reconstruction 



FIG. 16: Spatial tomography of the classical laser mode. 
Reflecting glass plate GP and cylindrical lenses CLl, CL2, 
create a one-dimensional field with two amplitude peaks in 
the 2 = plane. Cylindrical lens CL3 (oriented 90° from CLl 
and CL2) is varied in position, and intensity profiles in the 
z = D plane are imaged and recorded using spherical lens SL 
and the camera. Profiles are measured for 32 combinations 
of distances d and D. The piezoelectric transducer (PZT) 
introduces partial coherence between the two peaks. From 
McAlister et al. (1995). 



The non-interferometric method just described is best 
performed with an array detector to image the proba- 
bility distributions at different propagation distances^^. 
In addition, for a reliable reconstruction, it is necessary 
to ensure that the beam waist (its region of minimum 
spatial extent) occurs well within the measured zone. If 
this is not the case, then only partial state reconstruc- 
tion is possible. In the case of a limited scan, the density 
matrix in momentum representation can be measured ev- 
erywhere except for a band around the diagonal, whose 
width decreases as a larger range of longitudinal distances 
(time-of-flight) is measured (Raymer, 1997b). This was 
the case in a demonstration of transverse spatial QST 
of an ensemble of helium atoms (Kurtsiefer et al., 1997; 
Janicke and Wilkens, 1995). 

By using lenses (for light or for atoms) , the waist region 
can be brought into range for imaging, thus ensuring a 
reliable reconstruction (Raymer et al., 1994; McAlister 
et al, 1995). Suppose a beam propagates through a lens 
ai z = d, with focal length /, to a detection plane at 
z = D. The Maxwell- field propagator (73 1, in paraxial 
approximation, takes the form 



G{x' , x;t) — C exp 



ih{x) -f ifco 



L 



2Rc 



(76) 

where L = [D ~ d){l H- d/i?o), Rc ^ Ro + d, and R^^ = 
{D — d)^^ — f^^. Here C is a constant and h{x) is an 
unimportant phase function. 

The first measurements of this type were carried out 
for macroscopic ( "classical" ) fields from a laser (McAlis- 
ter et al., 1995). Figures 16 and 17 show the setup and 



the reconstructed data in the object plane, for the case 
of a two-peaked field distribution created by reflecting 
the signal beam from a two-sided reflecting glass plate. 
The field correlation function {E{xi)E* {X2)) (which in 
this case should not interpreted as a quantum density 
matrix) was reconstructed using the method in Raymer 
et al. (1994) — the inverse Radon transform of intensity 
distributions measured for different lens and detector po- 
sition combinations. 



For the data in Fig. 17 'a), the neutral-density filter ND 
was replaced by a beam block, so two coherent beams 
comprised the signal. Four lobes are seen in the recon- 
structed field correlation function. Figure 17 'c) shows 
the reconstructed field and the phase profile. For the 



Scanning a single detector would be prohibitive. 



data in Fig. [I7[b), the beam block BB was inserted in 
the lower beam (as shown) and the second beam compo- 
nent was created by reflection from a mirror M3 mounted 
on a translator driven by a random voltage, so two mutu- 
ally incoherent beams comprised the signal. In this case, 
the off-diagonal lobes are missing in the reconstructed 
field correlation function, as expected. This experiment 
verified the method of phase-space tomography for re- 
constructing spatial field correlations at the macroscopic 
level. Hansen (2000) applied a similar technique to recon- 
struct the optical mode emerging from a two-slit inter- 
ferometer and obtained a Wigner function with negative 
values, similar to that of Kurtsiefer et al. (1997). The 
method has been applied to study light scattering from 
complex fluids (Anhut et al., 2003). 

Other methods for classical wave-front reconstruction 
have since been developed (laconis and Walmsley, 1996; 
Cheng et al., 2000; Lee et al., 1999) and applied (Cheng 
and Raymer, 1999; Lee and Thomas, 2002; Reil and 
Thomas, 2005). A method has been elaborated for char- 
acterizing single-photon states in terms of a discrete spa- 
tial basis (Sasada and Okamoto, 2003; Langford et al.. 
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FIG. 17: Spatial tomography results for the classi- 
cal laser mode. Equal-separation contours showing the 
magnitude-squared of the reconstructed field correlation func- 
tion p(XiX')i foi' i^) the fully coherent field and (b) the par- 
tially coherent field, (c) The intensity profile (solid curve) and 
phase profile (dashed curve) of the reconstructed complex- 
wave field obtained in the fully coherent case. Axes are scaled 
transverse position x = ^/^o, where xq is a characteristic 
length. From McAhster et al, (1995). 
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FIG. 18: Top view of all-reflecting Sagnac interferometer. 
All mirrors are planar, including beam splitter BS. The beam 
travels out-of plane to reach the center mirror in the top- 
mirror configuration (surrounded by dashed lines), which is 
above the plane of the others. The displacement and tilt 
of the external steering mirror selects the phase-space point 
at which the Wigner distribution is measured. Signals from 
photon-counting detectors Dl and D2 are subtracted. 01 is 
an optical isolator for directing the reflected signal to D2. 



2004) and a proposal has been made for generalizing this 
to arbitrary beams (Dragoman, 2004). 



Interferometric reconstruction by wave-front 
inversion 



A technique for continuous-spatial-variable character- 
ization of single-photon fields was proposed by Mukamel 
et al. (2003), and recently implemented by Smith et. al. 
(2005). The method uses a parity-inverting Sagnac in- 
terferometer to measure the expectation value of the par- 
ity operator 11,^^ which, as first shown by Royer (1977), 
is proportional to the Wigner distribution at the phase 
space origin: 



W{Q, 0) = -Tr 



(77) 



The Wigner function at an arbitrary phase-space point 
can be determined by measuring the parity expectation 
value of the mode after the latter is displaced in the phase 
space in a manner similar to that proposed by Banaszek 
(1999) and discussed in the end of Sec. IIIB 5 for Wigner 
functions in the field quadrature space: 



1, 



-Tr 



(78) 



The parity operator's eigenstates are those with even and odd 
wavefunction, which correspond, respectively, to the eigenvalues 
1 and —1 



27 



Experimentally, the displacement D is implemented by 
physically shifting the mode location by a' and tilting its 
propagation direction by kx (see Fig. 18 1. 

The mode parity is measured as follows. One decom- 
poses the signal field into a sum of even and odd terms, 
E{x) — Ee{x) + Eo{x). Then the Wigner distribution 
( 78 1 evaluates to 



\E,{x')f) ~ (\E,{x')\ 



(fa 



W{x, kx) = 

(79) 

the terms in angular brackets being the experimentally 
measurable mean intensities or photon count rates for a 
given shift and tilt (x, kx)- This measurement is achieved 
by means of a dove-prism (Mukamel et at., 2003) or 
an all-reflecting (Smith et at., 2005) Sagnac interferom- 
eter as shown in Fig. [TSj The beam is split at beam 
splitter BS, after which the two beams travel in dif- 
ferent directions around the Sagnac loop. Each beam 
travels out-of plane to reach the center mirror in the 
top-mirror configuration, which has the effect of rotat- 
ing the wave fronts by ±90°, depending on direction, in 
the x-y plane. The net result is the interference of the 
original field with its (two-dimensional) parity-inverted 
image. Any odd-parity beam [E{-~x,—y) = —E{x.,y')\ 
passes through to detector Dl, while any even-parity 
beam \E{—x^ —y) = E{x, y)] reflects back toward the 
source, and is detected by D2. By subtracting the average 
count rates integrated over detector faces large enough to 
capture all signal light, one measures the Wigner distri- 
bution at a point in phase space, according to Eq. (79 1. 



One can also use only one detector, in which case the av- 
erage counting rate, as a function of x and kx, is propor- 
tional to the Wigner distribution plus a constant, which 
must be subtracted. 

In order to apply this technique in the photon- 
counting regime, one would like to use high quantum- 
efficiency avalanche photodiodes (APDs) operating in 
Geiger mode. Unfortunately, these typically have very 
small detector area (0.1 mm diameter), making them 
unsuitable for detecting beams with large intrinsic di- 
vergence. The experiment by Smith et al. (2005) used a 
single large-area photon-counting detector Dl. The de- 
tector was a photomultiplier tube with 5 mm diameter, 
11% efficient at wavelength 633 nm. 

Figure [19] shows results obtained for an expanded laser 
beam after passing through a single-sht or double-slit 
aperture placed in the beam just before the steering mir- 
ror. The beam was attenuated so that only a single pho- 
ton was typically present at any given time. Again, these 
results are similar to those obtained by Kurtsiefer et al. 
(1997) for a beam of helium atoms. Note that although 
the Wigner functions shown in Fig. [19] are the "Wigner 
functions of the single-photon Fock state", they repre- 
sent quantum objects fundamentally different from that 
plotted in Fig. [9] The latter describes the quantum state 
of a speciflc electromagnetic oscillator while the former 
describes the superposition of electromagnetic oscillators 
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FIG. 19: Measured Wigner functions for a slightly diverging 
beam passed through a single slit (left column) or double slit 
(right column) . Each shows a shear associated with beam di- 
vergence. In both cases, the interference fringes oscillate pos- 
itive and negative, as expected for a nonclassical momentum 
state. In the case of two slits, the fringes can be understood as 
resulting from a superposition ( "Schrodinger cat" -like) state 
of two well-separated components. From Smith et al. (2005). 



carrying a specific optical state. 

The ability to measure quantum states or wave func- 
tions for ensembles of single-photon states can be gener- 
alized to two-photon states. As pointed out in Mukamel 
et al. (2003) and Smith et al. (2005), if a photon pair in a 
position-entangled state pab can be separated, then each 
can be sent into a separate Sagnac interferometer, and 
subsequently detected. The rate of coincidence counts 
is proportional to a sum of terms, one of which is the 
two- photon Wigner distribution. 



W{xA,kxA,XB,kxi 



1 



TY[D-^pABDIlAflB], (80) 



[where D = D(xa, kxA, ^b, ^xb)] which can be extracted 
from the counting data. 

If the state is pure, this Wigner function can be trans- 
formed into the "two-photon wave function" , defined by 
(Smith and Raymer, 2007) 



^{xa,xb) = {y&c\E'^+\xa)E^+\xb)\^ 



AB) 



(81) 



Measuring the two-photon wave function would provide 
a complete characterization of position-entangled states, 
which are of interest in the context of Einstein- Podolsky- 
Rosen correlations or Bell's inequalities with photon po- 



^® A related quantity, the two-photon coincidence-detection ampli- 
tude, has been discussed in Nogueira et al. (2002), Walborn et 
al (2004), Keller and Rubin (1997), and Scully and Zubairy 
(1997). 
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sition and momentum variables (Howell et ai, 2004; 
Yarnall et ai, 2007). 



VII. SUMMARY AND OUTLOOK 

Prior to the beginning of the present century, quantum 
physics of light has been developing along two parallel 
avenues: "discrete-variable" and "continuous-variable" 
quantum optics. The continuous-variable community 
dealt primarily with the wave aspect of the electromag- 
netic field, studying quantum field noise, squeezing and 
quadrature entanglement. Homodyne detection was the 
primary tool for field characterization. The discrete- 
variable side of quantum optics concentrated on the par- 
ticle aspect of light: single photons, dual-rail qubits, and 
polarization entangled states. These objects were usually 
measured with single-photon detectors. 

These two aspects of quantum optics had little overlap 
with each other in terms of methodology, but experienced 
significant mutual infiuence. Novel results in the discrete- 
variable domain, such as demonstration of entanglement, 
quantum tomography, quantum teleportation, etc., were 
frequently followed by their continuous-variable analogs 
and vice versa. 

Theoretically, the difference between these two do- 
mains boils down to the choice of the basis in which states 
of an optical oscillator are represented: either quadrature 
(position or momentum) or energy eigenstates. From the 
experimental point of view, parametric down-conversion, 
the workhorse of quantum optical state production, can 
operate in either the weak or strong pumping modes. In 
the former case, we obtain discrete photon pairs, in the 
second, squeezing or quadrature entanglement. 

The division of quantum optics is thus caused not by 
fundamental but by pragmatic reasons. It is just that, 
until recently, our technology allowed us to generate only 
two classes of quantum states, giving us access to two 
small islands in the vast ocean of the optical Hilbert 
space. 

Developments of the last decade allowed us to over- 
come this separation. By applying a traditionally 
continuous-variable quantum characterization method 
(homodyne tomography) to discrete- variable quantum 
states (photons and qubits), researchers have constructed 
the bridge between the two islands, and then extended 
it by engineering and characterizing quantum states 
that belong to neither domain — such as displaced 
and photon-added states, squeezed Fock states and 
Schrodinger cats. 

In this review, we covered technological developments 
that led quantum optics to this breakthrough, placing a 
particular accent on continuous-wave tomography. We 
discussed new state-reconstruction algorithms, the tech- 
nology of time-domain homodyne detection, preparation 
of high-purity photons and qubits, and methods of quan- 
tum state engineering. We also reviewed methods of 
characterizing the modal structure of a quantum-optical 



state. 

Extrapolating the last years' results into the future, we 
can isolate certain open problems and future directions 
along which the field can be expected to develop. 

Reliable state-reconstruction algorithms. We have fo- 
cused on maximum likelihood estimation (MaxLik) be- 
cause it is straightforward to implement and offers im- 
provements over the inverse-linear-transform techniques 
such as inverse Radon. However, MaxLik is probably 
not the last word in QST algorithms. It can under- 
perform if only a small amount of data is available. 
In some cases this technique can yield zero probabili- 
ties for certain state components, which are not justified 
(Blume-Kohout, 2006). We predict that future recon- 
struction algorithms will combine maximum-likelihood 
with maximum-entropy and Bayesian methods (Fuchs 
and Schack, 2004). One attempt at such integration has 
already been reported (Rehacek and Hradil, 2004). 

Within MaxLik itself, particularly in application to 
OHT, there are a number of open questions. To what 
extent does a bias in the tomography scheme (non-unity 
sum of the POVM elements, see Sec. HI B 5 ) infiuence the 



reliability of state reconstruction? What is the optimal 
point for truncating the Hilbert space that would allow 
sufficiently complete but noise-free reconstruction? Does 
there exist a simple and reliable method for evaluating 
errors in quantum state estimation? 

Faster, low-noise homodyne detectors. As discussed 
in Sec. |IV A| there is a compromise between the band- 
width of the homodyne detector and its signal-to-noise 
ratio. Detectors with higher bandwidths can accom- 
modate higher laser repetition rates, permitting acqui- 
sition of larger data sets and eventually analysis of more 
complex states of light. Future study in this area is 
well-deserved, also given applications of faster homodyne 
detectors is continuous-variable quantum cryptography, 
with a promise of significant secret key transfer rate en- 
hancement. 

Applications of OHT in discrete quantum-information 
processing. Most of the optical protocols tested so far 
employed dual-rail qubits as quantum information car- 
riers (Kok et al., 2007). Accordingly, photon counting 
has been the method of choice for state measurement. 
As discussed in the Introduction, homodyne tomography 
provides much more complete information about a state 
of light (and thus performance of a quantum gate), but 
it is not yet commonly employed due to the relative com- 
plexity of its implementation. A goal for future research 
would be to simplify basic elements of homodyne de- 
tection — mode matching, local oscillator phase varia- 
tion, data acquisition — to an extent that would make 
OHT not much more complicated than photon counting. 
Homodyne tomography should then be applied for char- 
acterizing complex discrete-variable quantum states and 
protocols. Perhaps one of the first steps would be char- 
acterization of an entangled state of two dual-rail qubits 
— akin to that in James et al. (2001), but using OHT. 

Another important Ql-related application is testing 
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protocols on interfacing quantum information between 
light and stationary media. Examples of such are the 
experiments of Julsgaard et al. (2004) and Appel et 
al. (2008), which utilized homodyne tomography to study 
quantum memory for light. Homodyne measurements 
on light transmitted through an atomic ensemble per- 
mit tomographic reconstruction of its collective spin state 
(Feenholz et al., 2008), which is useful, for example, 
for characterizing quantum tomography from light onto 
atoms (Sherson et al., 2006). 

Continuous-variable process tomography. While mea- 
suring superoperators associated with a certain quantum 
process has been investigated theoretically (Chuang and 
Nielsen, 1997) and experimentally (Altepeter et al., 2003) 
for discrete variables for quite some time, the progress in 
the continuous- variable domain has been non-existent. 
This seems to be an important open problem, whose so- 
lution holds a promise to provide much more complete 
data about quantum processes than current methods. 

Quantum-state engineering, i.e. synthesis of arbitrary 
quantum states of light using nonclassical primitives 
(squeezed or Fock states), linear optics and conditional 
measurements. There exist a number of proposals for 
tackling this objective (reviewed in detail in Dell'Anno 
et al, 2006), for example by using coherent displace- 
ments and photon subtraction operations (Dakna et al., 
1999a,b; Fiurasek et al., 2005), repeated parametric 
down-conversion (Clausen et al., 2001) and continuous- 
variable postselection (Lance et al., 2006). To date, we 
have mastered quantum state engineering at the single- 
photon level: we can create any linear combination of the 
vacuum and single-photon Fock state. The next step is 
to bring this to the two-photon level. This can be done, 
for example, by applying modified photon addition oper- 
ations (Sec. VDl to single-rail qubits. 

Quantum optical engineering, as well as any other com- 
plex manipulation of light, requires high quality of the 
"raw material", i.e. initial squeezed and Fock states. 
Here we can see two possibilities for progress. On one 
hand, parametric down-conversion sources need to be 



improved to generate spectrally and spatially unentan- 
gled signal and idler photons, as well as pulsed squeez- 
ing in a single spectral mode (Wasilewski et al., 2006). 
On the other hand, it would be great to eliminate down- 
conversion altogether and employ solid-state, on-demand 
sources (Grangier et al., 2004). At present, such sources 
compromise between efficiency and spatiotemporal pu- 
rity, and thus cannot be employed in scalable quantum- 
optical engineering. We hope that the situation will 
change in the near future. Additionally, there may exist 
a possibility for improving the efficiency of such sources 
by means of linear optics and conditional measurements 
(Berry et al, 2006, 2007). 

Fundamental tests and new quantum protocols that are 
not restricted by either discrete- or continuous domains of 
quantum optics. Examples are loophole-free nonlocality 
tests (Garci'a-Patron et al., 2004; Nha and Carmichael, 
2004) and purification of continuous-variable entangle- 
ment (Opatrny et al, 2000; Browne et al, 2003). All 
"building blocks" of these protocols have already been 
experimentally demonstrated, but a task to put them to- 
gether in operational setups remains on the agenda. 

In summary, more work is needed before we gain full 
control over the optical Hilbert space. It is however worth 
the effort: if we have seen so many wonders within the 
boundaries of the two small islands colonized so far, who 
can predict what surprises await us in the vast expanses 
of the whole ocean? 
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